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SECTION  I 


INTRODUCTION 


Titanium  and  titanium  alloys  are  used  extensively  in  aerospace  systems  because  of  their 
high  strength-to-density  ratio,  temperature  resistance  and  excellent  corrosion  behavior. 
Titanium  alloys  have  been  widely  applied  to  aircraft  structural  components,  and  to  turbine 
engine  and  compressor  parts.  As  titanium  e'loys  have  found  application  in  different  environ- 
ments, especially  at  increasingly  higher  temperatures  and  pressures,  a renewed  interest  has 
developed  in  their  ignition  and  combustion  behavior  in  air. 

Within  the  last  ten  years,  the  incidence  of  titanium  fires  in  gas  turbine  engines  has 
markedly  increased.  Most  frequently,  catastrophic  massive-metal  fires  have  been  observed 
during  test  cell  performance  evaluations.  A number  of  factors  are  involved  in  this  trend;  most 
notably,  the  increase  in  thrust-to-weight  ratio  of  current  engines.  Lower  weight  demands  a more 
flexible  engine  and  lighter  components  and  these  factors,  when  added  to  the  desire  for  more 
thrust,  result  in  increased  pressure,  temperature  and  flow,  providing  an  environment  conducive 
to  ignition,  self-sustained  combustion  and  propagation  of  a titanium  fire. 

To  date,  the  primary  remedy  has  involved  substitution  of  steels  for  selected  titanium  parts, 
with  a resultant  weight  penalty,  and  increase  in  the  rotor/stator  clearances  with  a resulting 
performance  penalty.  As  these  penalties  may  be  unacceptable  in  certain  “performance-driven” 
military  systems,  the  Air  Force  has  initiated  a program  to  understand  the  phenomena  and 
develop  materials  approaches  to  eliminate  catastrophic  failure. 

The  combustion  cf  titanium  has  been  extensively  investigated,  with  the  vast  majority  of 
studies  concerned  with  powdered  metals,  wires  and  single  droplets.  Wnite  and  Ward  (Reference 
1)  provide  a comprehensive  literature  review  of  the  ignition  of  metals  in  oxygen  and  mixed-oxygen 
atmospheres  prior  to  1966.  The  most  notable  of  these  earlier  studies,  insofar  as  this  effort  is 
concerned,  is  the  work  of  Liftman,  et  al.,  (Reference  2)  which  is  one  of  the  few  studies  concerned 
with  the  spontaneous  ignition  of  massive  metal.  The  dependence  of  fracture  initiated  ignition  on 
pressure  and  atmosphere  reported  by  Liftman,  et  al.,  has  been  used  to  set  the  limits  for  titanium 
usage  in  gas  turbine  engines  as  recently  as  early  1976.  This  has  been  the  trend,  even  though 
fracture-initiation  may  be  the  least  important  source  of  ignition  in  the  modem  compressor 
environment. 

With  the  advent  of  high  energy  CO,  lasers,  several  investigators  have  noted  laser-initiated 
combustion  of  titanium  and  its  alloys.  Most  of  these  studies  have  investigated  the  initial 
conditions  and  dependence  of  laser  burnthrough  times  on  such  variables  as  beam  size,  Mach 
number,  velocity,  or  specimen  geometry.  More  recently,  the  Air  Force  has  initiated  several 
studies  to  consider  various  types  of  ignition,  such  as  rubbing  or  rapid  heating,  and  to  investigate 
and  model  self-sustained  combustion  of  titanium  alloys  in  the  gas  turbine  engine  compressor 
environment.  Wolf  (Reference  3),  has  investigated  the  isothermal  oxidation  kinetics  of  several 
titanium  alloys  and  the  conditions  which  lead  to  ignition  of  thin  sheets.  With  heating  rates  on  the 
order  of  100°C/sec  unalloyed  titanium  sheet  has  been  ignited  at  approximately  500°C  in  200  Torr 
stagnant  oxygen.  Clark,  et  al.  (Reference  4 and  5),  are  investigating  the  combustion  mechanism 
of  a variety  of  titanium  alloys,  as  part  of  a continuing  program  at  the  Bureau  of  Standards.  This 
excellent  work  is  providing  a first  definition  of  the  basic  metal  combustion  process  as  well  as 
temperature  measurements  during  combustion. 

To  date,  there  has  been  little  effort  to  understand  the  ignition  and  propagation  of  a titanium 
fire  on  massive  blade-like  shapes  in  flowing  gas.  Given  a rectangular  or  airfoil  shaped  titanium 
alloy  in  a compressor  environment  and  sufficient  energy  for  ignition  at  a corner  or  edge,  it  is 
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expected  that  several  factors  will  determine  whether  propagation  occurs.  The  combination  of  gas- 
temperature,  pressure  and  velocity  control  the  concentration  of  reactants  at  a combusting 
surface,  the  aerodynamic  shear  for  removal  of  melted  materials,  and  the  convective  cooling  of  the 
surface.  No  prior  efforts  have  determined  the  dependence  of  ignition  or  propagation  on  these 
factors.  The  objective  of  this  effort  was  to  develop  an  understanding  of  the  interrelationships  of 
the  important  factors  governing  ignition  and  combustion  propagation  over  a broad  range  of 
environmental  conditions. 

The  efforts  described  in  this  report  include  an  experimental  evaluation  of  the  geometrical 
and  environmental  effects  on  the  combustion  of  Ti  8-1-1  airfoils,  and  the  development  of  an 
analytical  model  for  the  prediction  of  the  ignition  and  subsequent  combustion  of  titanium  and  its 
alloys. 


SECTION  II 


PROGRAM  SUMMARY 


A comprehensive  program  has  been  conducted  to  develop  an  in-depth  understanding  of  the 
interrelationship  of  the  factors  governing  the  ignition-propagation  sequence  of  titanium  over  a 
broad  range  of  environmental  conditions,  including  those  of  current  and  advanced  gas  turbine 
engine  operations.  Included  in  the  program  were  an  experimental  investigations  of  titanium 
combustion,  and  an  extensive  analytical  effort  to  develop  a predictive  model  for  evaluating  the 
susceptibility  of  titanium  alloys  to  ignition  and  subsequent  propagation. 

In  the  experimental  phase  of  the  program,  73  combustion  tests  were  performed  with 
titanium  8-1-1  specimen  airfoils  over  a wide  range  of  temperature,  pressure,  and  airstream 
velocity.  Combustion  was  initiated  with  a laser  ignition  source,  and  the  ignition-combustion 
sequence  for  each  test  was  recorded  with  high  speed  cinematography.  Subsequent 
phenomenological  analysis,  and  frame  by  frame  differential  motion  analysis  provided  the  basis 
for  development  of  the  analytical  models. 

In  the  initial  analytical  phase,  a phenomenological  sequence  was  described,  based  on  the 
observed  test  results,  and  a model  was  proposed  to  explain  the  observed  phenomena.  The  model 
was  subsequently  developed,  programmed  for  computer  analysis,  and  compared  with  test  data  to 
evaluate  its  validity.  The  model  appeared  to  be  qualitatively  in  agreement  with  test  data,  but 
quantitative  disagreement  of  certain  aspects  of  the  model  indicated  the  need  for  further 
refinement.  In  particular,  liquid  melt  spread-rates  over  the  airfoil  and  spanwise  spreading  of  the 
burning  zone  were  inconsistent  with  observed  data. 

In  the  advanced  development  phase  of  the  analytical  effort,  the  model  was  generalized  to 
include  centrifugal  force  effects,  thus  allowing  analysis  of  rotating  compressor  blades.  A detailed 
model  was  developed  for  the  flow  of  melt  and  for  the  transfer  of  heat  and  oxygen  within  the 
flowing  melt.  In  addition,  physical  property  models  have  been  modified  to  incorporate  the  effect 
of  dissolved  oxygen,  and  an  aerodynamic  loss  model  has  been  incorporated  in  the  liquid  flow 
equations  to  account  for  stripping  of  melt  from  the  burning  surface. 
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EXPERIMENTAL  PROGRAM 


An  experimental  program  was  conducted  to  provide  supporting  data  for  development  of  the 
analytical  model.  A total  of  73  tests  were  conducted,  of  which  44  resulted  in  burn  damage  to  the 
airfoil  specimen. 

A.  TEST  FACILITIES 

The  titanium  combustion  test  rig  used  for  these  studies  is  a small  wind  tunnel  driven  from 
a 241  MPa  (350  psi)  compressed  air  supply.  Ancillary  conditioning  equipment  permits  the 
simulation  of  a wide  range  of  environmental  combinations  of  air  pressure  — up  to  0.97  MPa  (140 
psia),  temperature  — up  to  550°C  (1022°F)  and  velocity  — up  to  335  m/sec  (1100  ft/sec).  Test 
specimens  are  mounted  in  a rectangular  test  chamber.  Instrumentation  is  provided  to  monitor 
pressure  and  temperature  at  strategic  locations  in/on  the  rig.  The  overall  arrangement  of  the  test 
rig  is  shown  in  Figure  1. 


Pq  = Orifice  Downstream  Pressure 
Pq  =Test  Chamber  Pressure 
Tq  = Orifice  Air  Temperature 
Ts  = Rig  Skin  Temperature 
Tq  =Test  Chamber  Temperature 
Tow  = Window  Temperature 

FT)  139R6* 

Figure  ].  Titanium  Combustion  Test  Rig 

Air,  supplied  from  a large  compressor,  is  passed  through  a gas-fired  indirect  heater  and  a 
flow-measuring  orifice  prior  to  entering  the  test  section.  The  test  section,  shown  schematically  in 
Figure  2 is  a 19.1  mm  by  50.1  mm  (0.75  in.  by  2.0  in.)  rectangular  channel  with  a bellmouth  inlet 
and  75  mm  (3.0  in.)  of  straight  section  upstream  of  the  airfoil  leading  edge.  The  test  specimen 
mounts  in  a cylindrical  carrier  which  is  inserted  into  the  test  section.  (See  Figure  3.)  Orifice 
plates,  upstream  and  downstream  of  the  test  section,  provide  control  of  the  flowrate  and  pressure 
level.  Thermocouples  are  positioned  to  provide  temperatures  at  the  flow-measuring  orifice,  in  the 
specimen  test  chamber,  on  the  rig  skin,  and  at  the  laser  and  camera  windows.  Points  are  provided 
to  measure  flow-control  orifice  upstream  and  downstream  pressure  as  well  as  static  pressure  in 
the  test  chamber.  (See  Figure  1.)  Airs tream  velocity  is  determined  by  calculation  using  the 
differential  pressure  across  the  flow-control  orifice. 
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Figure  3.  Test  Specimen  in  Holder 

The  test  section  has  two  windows,  one  for  laser  irradiation  of  the  specimen  and  one  for 
photographic  and  video  viewing  of  the  test  specimen . These  windows  are  mounted  in  a port  on  the 
side  of  the  rig  approximately  20.3  cm  (8  in.)  from  the  test  specimen.  The  window  for  the  camera 
is  optically  flat,  fused  quartz  which  is  63.5  mm  (2.500  in.)  in  diameter  and  12.7  mm  (0.500  in.) 
thick.  The  laser  beam  window  is  38.1  mm  (1.500  in.)  in  diameter  and  6.35  (0.250  in.)  thick  zinc 
selenide.  This  material  offers  excellent  transmissivity  for  the  10.6-micron  wavelength  emission  of 
the  CO,  laser  beam.  The  zinc  selenide  window  has  an  antireflective  coating  on  both  faces  to 
minimize  reflection  and  scatter  of  the  beam.  The  quartz  window  has  excellent  optical  clarity  for 
visible  light  but  will  not  transmit  the  CO,  laser  light,  thus  protecting  the  camera  lens  from 
scattered  reflections  of  the  laser  beam.  The  windows  are  protected  from  the  high-temperature  test 
environment  by  a water  jacket  to  absorb  conducted  heat  in  the  metai  housing  and  an  air  injection 
system  to  film-cool  the  optical  surfaces. 
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The  arrangement  of  the  test  section  and  other  supporting  test  equipment  is  shown  in  Figures 
4 and  5.  Because  of  environmental  requirements,  the  laser  equipment  is  located  in  the  air- 
conditioned  control  room  and  the  beam  passed  through  a port  in  the  concrete  blast  wall.  A high- 
speed Hycam  motion  picture  camera  and  a video  camera  with  a tape  recording/playback  system 
are  arranged,  through  the  use  of  a beamsplitter,  to  permit  simultaneous  photographic  recording 
and  real-time  video  observation.  The  energy  required  for  specimen  ignition  is  supplied  by  a CRL 
Model  41  laser.  This  electric -discharge,  water-cooled  CO,  laser  system  is  capable  of  providing  an 
output  of  approximately  250  watts  in  the  TEMoo  mode  at  a transmission  frequency  of 
10.6  microns. 

The  beam  is  defocused  at  the  specimen  to  a diameter  of  approximately  2 mm  to  yeild  an 
incident  average  power  density,  after  transmission  losses,  of  approximately  2.5  Kw/cm1  absorbed 
by  the  specimen.  A coincident  helium-neon  laser  is  used  to  provide  a visible  red  beam  for 
alignment  of  the  CO,  laser  beam  on  the  titanium  specimen.  This  red  alignment  beam  is  observed 
using  the  videotape  system. 

B.  TEST  RUN  PROCEDURE 

Prior  to  a test  run,  the  test  specimen  is  secured  in  the  specimen  holder  and  mounted  in  the 
test  chamber.  Run  values  for  air  velocity  and  chamber  pressure  are  established  by  installing  a 
specific  combination  of  flow-control  and  backpressure  orifices  in  the  test  rig.  After  setting  the 
heater  to  the  desired  temperature,  airflow  is  introduced  into  the  system  and  the  test  chamber 
allowed  to  come  to  thermal  equilibrium.  Near  thermal  equilibrium,  the  laser  optics  are  aligned 
using  the  visible  helium-neon  laser  beam  and  the  videotape  system.  Specimen  illumination 
during  laser,  video  and  high-speed  camera  alignment  and  focusing  is  provided  by  a high  intensity 
light  transmitted  from  its  source  by  fiber  optics.  During  the  run,  light  from  the  specimen  ignition 
and  burning  is  sufficiently  intense  to  permit  photodocumentation  of  burn  propagation  and  melt 
transportation.  An  event  marker  (light  pulse)  is  recorded  on  the  side  of  the  film  on  all  runs  to 
annotate  the  start  and  finish  of  the  laser  action  during  the  run.  The  high-speed  films  are  also 
marked  with  light  pulses  from  a 1000-Hz  timing  generator  to  provide  an  absolute  time  reference 
for  event  sequences. 

Just  prior  to  the  start  of  a run,  final  temperature  and  pressure  adjustments  are  made  by 
judicious  throttling  of  the  hot  air  valve  and  a cold  air  bleed  input  valve.  Simultaneously,  a final 
adjustment,  if  required,  is  made  to  the  laser  alignment.  Final  run  parameters  are  then  recorded 
and  the  test  is  initiated  by  a time-sequenced  switch  which,  when  actuated,  starts  the  high-speed 
camera.  Approximately  2 sec  after  camera  start  the  sequencer  opens  the  laser  shutter  to  irradiate 
the  specimen,  thereby  starting  the  run.  The  laser  remains  on  for  5 sec  before  the  sequencer  closes 
the  shutter.  This  time  can  be  manually  overridden  when  ignition  occurs  before  the  5 sec  has 
elapsed.  The  camera  is  allowed  to  run  out  of  film  (about  16  sec)  and  is  sequenced  off  at 
approximately  20  sec.  The  video  system  remains  on  at  all  times  in  a CCTV  mode.  Video  tape 
recording  is  controlled  manually  during  a run  sequence.  The  instant  replay  and  a slow- 
motion/stop-frame  capability  of  the  color  video  recording  permits  immediate  review  of  the  test 
run  for  its  potential  impact  on  the  next  test  run. 

C.  TEST  SPECIMENS 

All  combustion  tests  were  performed  with  a single  titanium  alloy.  Because  of  its  extensive 
use  in  gas  turbine  engine  applications,  Ti-8Al-lMo-lV  was  selected  to  provide  a baseline 
performance  upon  which  to  develop  the  analytical  model. 

The  basic  specimen  configuration  was  an  uncambered  airfoil,  aerodvnamically  represent- 
ative of  current  compressor  airfoils.  The  design  used,  shown  in  Figure  6,  is  an  uncambered  airfoil 


Concrete  Blast  Wall 


Arrangement  of  I. riser  Ignition  and  Photographic  Recording  System 


Test  Chamber 


Titanium  Combustion  Test  Rip  Setup 


I 


with  cylindrical  surfaces  of  115.6  mm  (4.55  in.)  radius,  rounded  leading  and  trailing  edges  with 
nominal  radius  of  0.076  mm  (0.003  in.),  and  chord  of  25.4  mm  (1.00  in.).  The  maximum  thickness 
is  approximately  1.27  mm  (0.050  in  ),  and  the  total  contoured  span  is  56.39  mm  (2.220  in.).  The 
specimens  were  machined  from  sheet  by  EDM,  polished  to  reduce  surface  roughness,  and  vapor- 
blasted  to  yield  a uniform  matte  surface.  All  specimens  were  then  profiled  to  determine  their 
exact  contour,  and  profiles  were  recorded  for  future  reference. 

Twenty -three  tests  were  performed  with  specimen  airfoils  in  the  vapor-blasted  condition. 
Because  of  apparent  difficulty  in  attaining  ignition,  the  remaining  specimens  were  coated  with  a 
black  nickel  coating,  increasing  absorptivity  to  near  6091,  thereby  maximizing  coupling  of  laser 
energy  to  the  specimens.  The  thin  nickel  plating  was  applied  only  on  a small  area  at  the  leading 
edge  ignition  site,  thereby  having  no  effect  on  subsequent  combustion.  Four  specimens  were 
prepared  as  flat  plates,  with  a knife-edge  on  one  side  to  facilitate  ignition,  as  shown  in  Figure  7. 

D.  EXPERIMENTAL  RESULTS 

A total  of  73  laser-initiated  combustion  tests  have  been  conducted  in  the  test  program,  of 
which  44  sustained  some  level  of  damage,  and  27  underwent  self-sustained  combustion.  No 
attempt  was  made  to  study  the  laser-initiated  ignition  process  in  detail,  the  laser  being  used  only 
to  provide  a reproducible  and  remote  ignition  source. 

Data  resulting  from  the  test  program  are  tabulated  in  Table  I.  Test  conditions  of  static 
pressure,  total  temperature,  free -stream  air  velocity  and  Mach  number  are  tabulated  with  results 
listed  as  material  combusted  or  lost  (‘t  damage).  Damage  is  defined  as  the  percent  of  exposed 
specimen  lost  (burned)  during  the  test,  excluding  the  portion  held  in  the  support.  As  combustion 
nears  the  support  (holder),  the  relatively  massive  metal  interferes  with  the  flow,  and  acts  as  a 
heat  sink,  always  resulting  in  a significantly  diminished  combustion  rate,  or  in  extinction. 
Consequently,  specimen  damage  above  approximately  80'’)  should  be  considered  total  self- 
sustained  combustion  of  the  metal.  No  attempt  was  made  to  account  for  this  behavior  by 
normalization  of  the  data. 

Figure  8 shows  a movie  frame  sequence  during  the  “steady  state”  portion  of  a typical  test. 
This  sequence,  beginning  after  40°r  of  the  chord  has  burned,  is  a flat  plate  specimen  being  tested 
under  a complementary  test  program  (under  Contract  N62269-77-C-0481).  The  bright  areas 
represent  molten  titanium  being  spread  over  the  airfoil  by  airstream  friction.  The  flowing  melt  is 
seen  to  spread  over  the  surface  much  more  rapidly  than  the  leading  edge  bums.  Occasional 
flashes  of  sparks  are  seen  in  the  background  as  melt  is  stripped  from  the  burning  surface  by  the 
airstream. 

Figure  9 shows  a time  sequence  composite  for  a typical  test,  developed  from  the 
cinematographic  data.  The  burning  leading  edge  is  seen  to  rapidly  progress  across  the  specimen 
in  the  chordwise  direction,  then  slowly  burn  in  the  spanwise  direction. 

The  effect  of  aerodynamic  angle  of  attack  on  the  rate  and  nature  of  airfoil  combustion  is 
shown  in  Figure  10.  Time  sequence  plots  of  the  burning  leading  edge  for  two  tests,  PW66  and 
PW67,  performed  at  the  same  environmental  conditions  but  with  different  airfoil  orientation  are 
shown  in  Figures  10  (a)  and  10  (b),  respectively.  With  parallel  flow  (zero  angle  of  attack), 
combustion  is  primarily  in  the  chordwise  direction,  typical  of  most  of  the  test  results.  The  effect 
of  inclining  the  airfoil  10  deg  relative  to  the  air  stream  direction  is  shown  in  Figure  10  (b).  The 
effect  of  the  difference  in  aerodynamics  on  combustion  propagation  is  readily  seen,  and  appears 
to  be  due  to  boundary  layer  separation  on  one  side  of  the  airfoil.  Observation  of  the  high  speed 
movies  of  Test  PW67  indicate  that  nearly  all  of  the  liquid  metal  flow  is  confined  to  the  airfoil 
pressure  surface.  Chordwise  propagation  appears  to  be  somewhat  faster  at  zero  angle  of  attack, 
but  total  propagation  and  the  amount  of  resulting  damage  appear  significantly  greater  for  the 
10  deg  angle  test. 
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'igure  6.  Titanium  Specimen  Configuration 
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Figure  7.  Flat  Plate  Specimen 


TABLE  I.  TEST  RESULTS 


Test  Condition Results 

P, T,  V'  Damage 

Run  No.  MPa  (ptia)  °C  (°F)  m/sec  (ft/sec)  Mn  % Loss 


Comments 


PW-1 

0.45 

(66) 

204 

(400) 

186 

(609) 

0.432 

0 

PW-2 

0.69 

(100) 

204 

(400) 

186 

(609) 

0.432 

0 

PW-3 

0.97 

(140) 

204 

(400) 

186 

(611) 

0.433 

0 

PW-4 

0.69 

(100) 

316 

(600) 

212 

(696) 

0.444 

0 

PW-5 

0.69 

(100) 

316 

(600) 

209 

(687) 

0.438 

0 

PW-6 

0.69 

(100) 

204 

(400) 

188 

(616) 

0.436 

0 

PW-7 

0.97 

(140) 

204 

(400) 

116 

(379) 

0.266 

0 

PW-8 

0.55 

(80) 

204 

(400) 

116 

(378) 

0.265 

0 

PW-9 

0.21 

(30) 

316 

(600) 

103 

(338) 

0.371 

0 

PW-10 

0.31 

(46) 

316 

(600) 

130 

(426) 

0.269 

0 

PW-11 

0.21 

(30) 

407 

(765) 

144 

(472) 

0.277 

0 

PW-12 

0.21 

(30) 

316 

(600) 

128 

(419) 

0.264 

0 

PW-13 

0.21 

(30) 

427 

(800) 

139 

(467) 

0.264 

0 

PW-14 

0.21 

(30) 

318 

(606) 

69 

(192) 

0.120 

0 

PW-15 

0.21 

(30) 

316 

(600) 

41 

(136) 

0.085 

0.2 

PW-16 

0.45 

(66) 

427 

(800) 

326 

(1066) 

0.637 

19.3 

PW-1 7 

0.45 

(66) 

316 

(600) 

306 

(1001) 

0.663 

0 

PW-18 

0.46 

(66) 

371 

(700) 

306 

(1000) 

0.649 

0 

PW-19 

0.45 

(66) 

427 

(800) 

324 

(1062) 

0.634 

0 

PW-20 

0.52 

(75) 

482 

(900) 

317 

(1039) 

0.595 

19.2 

PW-21 

0.52 

(76) 

467 

(865) 

313 

(1027) 

0.598 

15.6 

PW-22 

0.62 

(90) 

482 

(900) 

329 

(1079) 

0.619 

0 

PW-23 

0.86 

(125) 

316 

(600) 

276 

(906) 

0.586 

0 

PW-24 

0.62 

(90) 

427 

(800) 

304 

(998) 

0.593 

14.4 

PW-25 

0.62 

(90) 

371 

(700) 

290 

(961) 

0.589 

13.0 

PW-26 

0.62 

(90) 

371 

(700) 

291 

(966) 

0.592 

13.2 

PW-27 

0.61 

(89) 

371 

(700) 

286 

(935) 

0.579 

0 
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TABLE  I.  TEST  RESULTS  (Concluded) 


Teat  Condition Results 

P,  T,  V Damage 

Run  No.  MPa  (psia)  °C  (°F)  m/sec  (ft/sec)  Mn  % Loss Comments 


PW-28 

0.69 

(100) 

316 

PW-29 

0.45 

(66) 

316 

PW-30 

0.62 

(90) 

204 

PW-31 

0.45 

(66) 

204 

PW-32 

0.62 

(90) 

260 

PW-33 

0.62 

(90) 

316 

PW-34 

0.62 

(90) 

316 

PW-35 

0.45 

(65) 

371 

PW-36 

0.62 

(90) 

427 

PW-37 

0.45 

(65) 

260 

PW-38 

0.62 

(90) 

371 

PW-39 

0.21 

(30) 

371 

PW-40 

0.62 

(90) 

316 

PW-41 

0.62 

(90) 

427 

PW-42 

0.45 

(65) 

316 

PW-43 

0.62 

(90) 

316 

PW-44 

0.62 

(90) 

260 

PW-45 

0.45 

(65) 

260 

PW-46 

0.62 

(90) 

204 

PW-47 

0.62 

(90) 

149 

PW-48 

0.31 

(45) 

316 

PW-49 

0.21 

(30) 

371 

PW-50 

0.62 

(90) 

316 

PW-51 

0.45 

(66) 

371 

PW-52 

0.48 

(70) 

427 

PW-53 

0.97 

(140) 

316 

PW-54 

0.69 

(100) 

316 

PW-55 

0.62 

(90) 

427 

PW-56 

0.69 

(100) 

316 

PW-57 

0.97 

(140) 

204 

PW-58 

0.69 

(100) 

204 

PW-59 

0.45 

(65) 

427 

PW-60 

0.62 

(90) 

427 

PW-61 

0.62 

(90) 

427 

PW-62 

0.62 

(90) 

427 

PW-63 

0.62 

(90) 

427 

PW-64 

0.97 

(140) 

482 

PW-65 

0.97 

(140) 

316 

PW-66 

0.97 

(140) 

427 

PW-67 

0.97 

(140) 

427 

PW-68 

0.69 

(100) 

316 

PW-69 

0.69 

(100) 

316 

PW-70 

0.69 

(100) 

427 

PW-71 

0.69 

(100) 

427 

PW-72 

0.62 

(90) 

427 

PW-73 

0.62 

(90) 

371 

(600) 

278 

(911) 

0.590 

(600) 

278 

(912) 

0.591 

(400) 

248 

(813) 

0.585 

(400) 

247 

(811) 

0.583 

(500) 

265 

(871) 

0.593 

(600) 

279 

(915) 

0.593 

(600) 

248 

(815) 

0.525 

(700) 

290 

(950) 

0.589 

(800) 

302 

(991) 

0.589 

(500) 

262 

(858) 

0.584 

(700) 

261 

(856) 

0.526 

(700) 

290 

(951) 

0.589 

(600) 

126 

(412) 

0.260 

(800) 

137 

(449) 

0.260 

(600) 

123 

(405) 

0.255 

(600) 

121 

(398) 

0.251 

(500) 

119 

(392) 

0.260 

(500) 

117 

(385) 

0.255 

(400) 

113 

(371) 

0.260 

(300) 

106 

(349) 

0.260 

(600) 

121 

(398) 

0.251 

(700) 

118 

(386) 

0.232 

(600) 

208 

(661) 

0.435 

(700) 

216 

(708) 

0.432 

(800) 

226 

(740) 

0.433 

(600) 

210 

(688) 

0.439 

(600) 

209 

(687) 

0.439 

(800) 

226 

(743) 

0.435 

(600) 

206 

(677) 

0.432 

(400) 

ia5 

(607) 

0.430 

(400) 

185 

(608) 

0.430 

(800) 

225 

(738) 

0.430 

(800) 

226 

(743) 

0.435 

(800) 

226 

(743) 

0.435 

(800) 

226 

(743) 

0.435 

(800) 

226 

(743) 

0.435 

(900) 

44 

(143) 

0.079 

(600) 

38 

(126) 

0.079 

(800) 

46 

(150) 

0.086 

(800) 

46 

(150) 

0.086 

(600) 

199 

(652) 

0.415 

(600) 

206 

(677) 

0.432 

(800) 

244 

(799) 

0.469 

(800) 

244 

(799) 

0.469 

(800) 

271 

(889) 

0.525 

(700) 

288 

(945) 

0.585 

0 

10.7 
0 

0.2 

0.4 

12.2 

0.2  10  deg  Attack  Angle 

13.4 

13.8 
0.1 

0 10  deg  Attack  Angle 

0.4 

16.1 

21.2 

17.0 

35.7  10  deg  Attack  Angle 

9.3 
0.2 

6.0 
0.1 
0.1 
0 

0.5 

0.5 

16.3 
0 

0 

45.9 

13.7 

10.4 
0.3 

20.8 
0 

0.7 

0 

0 

,56.2 

1.0 

30.3 

82.6  10  deg  Attack  Angle 

1.1  10  deg  Attack  Angle 

10.4  0.97  mm  thickness 

64.5 

82.8  1.55  mm  thickness 

0.1  10  deg  Attack  Angle 

1.3 
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Figure  8.  Combustion  Propagation  Across  Specimen  Airfoil 


Figures  11  through  14  illustrate  the  effect  of  test  condition  on  the  chordwise  burning  rate. 
Burning  rates  are  observed  to  vary  strongly  with  time,  probably  indicating  a strong  dependence 
on  chordwise  location  and  amount  of  downstream  surface  (toward  trailing  edge)  covered  with 
melt.  For  direct  comparison,  data  should  be  considered  at  times  corresponding  to  equal  amounts 
of  melt  mass,  which  is  equivalent  to  identical  chordwise  location.  A composite  of  these  data  is 
shown  in  Figure  15,  with  burned  mass  (location)  as  a parameter.  For  the  limited  range  of  pressure 
and  temrerature  considered  in  this  figure,  burning  rate  is  observed  to  be  strongly  dependant  on 
airstream  velocity. 

The  effect  of  the  environmental  conditions  on  the  threshold  of  self-sustained  combustion  is 
illustrated  in  Figure  16.  In  this  figure  the  effect  of  pressure,  velocity,  temperature,  and  geometry 
is  combined  in  the  dimensionless  Reynolds  number  at  the  leading  edge,  based  on  the  leading  edge 
thickness  and  all  properties  evaluated  at  the  “melting  film”  temperature  (average  of  airstream 
and  titanium  melting  temperatures).  Allowing  for  a moderate  amount  of  data  scatter,  the  curve 
in  Figure  16  reasonably  differentiates  the  combustion  threshold.  The  illustrated  curve  is  based  on 
an  earlier  one-dimensional  analysis  of  the  combustion  process,  performed  at  P&WA  prior  to  the 
current  study. 
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Figure  1 1.  Leading  Edge  Burning  Rate , Test  PW  40 
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Figure  12.  Leading  Edge  Burning  Rate,  Test  PW  66 
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Figure  13.  Leading  Edge  Burning  Rate,  Test  PW  67 
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ure  14.  Leading  Edge  Burning  Rate,  Test  PW  71 
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Figure  15.  Relation  of  Burning  Rate  to  Airstream  Velocity 
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Figure  16.  Threshold  of  Self-Sustained  Combustion 


SECTION  IV 


BASIC  ANALYTICAL  MODEL 


To  provide  a procedure  for  predicting  the  fire  susceptibility  of  titanium  structures  in 
aerodynamic  environments,  an  analytical  model  has  been  developed.  Based  on  experimental 
observations,  a scenario  was  developed  for  the  ignition-combustion  process,  a physical  model  was 
assumed  to  represent  the  scenario,  and  an  analytical  model  was  formulated  to  characterize  the 
physical  model.  This  initial  (basic)  model  was  developed  for  static  systems,  with  no  consideration 
of  rotational  motion. 

A.  IGNITION  — PROPAGATION  SCENARIO 

In  order  to  clarify  the  relationship  of  combustion  rates,  total  damage  and  ignition,  it  is 
important  to  consider  the  entire  ignition-propagation  sequence  with  regard  to  the  local 
environment.  Based  on  observation  of  high  speed  movies  and  applied  to  a gas  turbine  engine 
environment,  the  following  scenario  can  be  developed  to  describe  the  titanium  fire  phenomena: 
A titanium  specimen  is  exposed  to  and  is  in  thermal  equilibrium  with  an  aerodynamic 
environment.  Heat  generated  by  normal  oxidation  is  balanced  by  losses  to  the  surrounding 
environment.  An  ignition  initiation  event  creates  an  imbalance  in  this  equilibrium  system.  Such 
an  event  may  take  one  of  several  forms,  i.e.,  heat  input  due  to  mechanical  friction,  thermal 
radiation,  aerodynamic  heating  or  fracture  which  produces  an  atomically  fresh  surface.  In 
response  to  one  of  these  events,  metal  temperature  increases  in  the  region  of  heat  input,  resulting 
in  increased  heat  generation  from  oxidation.  Heat  loss  to  the  environment  and  to  the  sample 
increases  accordingly  to  establish  a new  equilibrium  state.  If  thermal  equilibrium  is  reached  with 
the  local  metal  surface  temperature  below  its  melting  point,  ignition  does  not  occur,  and  the 
system  returns  to  its  original  state  when  the  ignition  source  is  removed.  If  melting  occurs,  the 
material  ignites  and,  with  removal  of  the  ignition  source,  may  cease  burning  if  the  thermal 
balance  is  insufficient  to  supply  the  necessary  latent  heat  to  continue  local  melting.  If  the  rate  of 
local  heat  generation  is  greater  than  the  rate  of  losses,  self-sustained  combustion  occurs.  In  the 
case  of  self-sustained  combustion  and  subsequent  propagation,  test  results  demonstrate  that  a 
significant  amount  of  molten  material  is  generated  at  the  leading  edge  of  the  specimen  and 
adheres  to  and  flows  over  the  surface  in  response  to  aerodynamic  forces.  This  flowing  liquid  with 
considerable  dissolved  oxygen  continues  to  react  as  it  flows,  thereby  maintaining  its  liquid  state 
and  providing  additional  energy  to  the  specimen  by  conduction.  The  adherence  of  this  molten 
metal  in  certain  environments  is  critical  to  the  maintenance  of  self-sustained  combustion,  and  is 
the  principal  mechanism  for  propagation  of  titanium  combustion  at  moderate-to-great  airstream 
velocities.  In  general,  the  transport  of  energy  along  the  surface  in  the  molten  flow  far  exceeds  the 
capability  for  heat  loss  by  conduction  with  the  solid  substrate.  A more  accurate  description  of  this 
process  must  come  from  analytical  modeling. 

B.  PHYSICAL  MODEL 

The  physical  model,  shown  in  Figure  17,  describes  the  essential  features  of  the  combustion 
of  a titanium  airfoil  in  an  airstream.  This  model  formulated  on  the  basis  of  photographic  evidence 
from  the  test  program,  assumes  that  molten  metal  originates  at  the  leading  edge  and  flows  on  the 
surface  in  the  direction  of  the  airstream,  driven  along  by  aerodynamic  shear  and  drag  forces. 


Part  of  the  liquid  is  carried  away  by  the  airstream,  the  remainder  being  retained  on  the 
surface,  where  it  continues  to  react.  Some  of  the  liquid  solidifies  as  it  flows  over  the  surface, 
yielding  its  latent  heat  to  the  airfoil.  As  the  process  continues,  the  solidified  “recast”  layer 
subsequently  melts  and  reacts,  along  with  the  original  airfoil  structure.  Heat  supplied  to  the 
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surface  from  external  sources,  plus  heat  flow  from  the  reacting  surface  to  the  underlying  solid 
substrate,  provides  the  energy  necessary  to  cause  propagation  of  the  process.  Simultaneously, 
heat  is  lost  to  the  surrounding  environment  by  convection  and  radiation,  retarding  the  process. 
Thermal  balance  of  the  heat  gains  and  losses  established  whether  the  burning  process  will  be  self- 
sustained,  and  at  what  rate  propagation  will  occur. 


Figure  17.  Physical  Model  of  Airfoil  Combustion 


C.  ANALYTICAL  APPROACH 


The  physical  model  depicted  in  Figure  17  has  been  analytically  simulated  using  a 
combination  of  lumped  mass  and  differential  numerical  techniques.  A lumped  mass  representa- 
tion of  a burning  airfoil,  with  the  various  heat  transfer  processes  indicated  is  shown  in  Figure  18. 
Each  lumped  mass,  or  node,  represents  a finite  fraction  of  the  airfoil,  and  is  presumed  to  be  in 
thermal  communication  with  its  surrounding  nodes  and  with  the  gas  flow  on  the  surface.  In 
addition,  if  liquid  metal  covers  the  surface  at  the  node  location,  the  temperature  of  the  liquid  is 
assumed  to  be  independent  of  the  temperature  of  the  nodal  mass,  but  heat  can  flow  between  the 
liquid  metal  and  the  nodal  mass. 


An  energy  equation  can  be  written  for  a nodal  mass  in  a generalized  form,  equating  the  sum 
of  all  sources  of  heat  to  the  rate  of  internal  energy  change.  The  energy  balance  takes  the  form: 


dE 

dt 


Qc.»nductlon  qronvpctlon  qrndlatlon 


(1) 


blnput  Qreacllon  **, 

The  rate  of  internal  energy  change,  dE/dt,  refers  to  the  energy  of  the  node,  and  energy 
change  can  be  in  the  form  of  sensible  heating  (temperature  change)  or  phase  change 
(melting/solidification).  Two  different  cases  are  noted:  (a)  no  liquid  on  the  surface,  and  (b)  the 
surface  covered  with  liquid  metal. 
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Figure  18.  Lumped  Mass  Thermal  Model 


For  the  first  case,  the  nodal  mass  represents  the  entire  cross-sectional  mass  of  the  airfoil. 
Conduction  occurs  to  neighboring  nodes,  convection  and  radiation  to  the  surrounding 
environment,  and  heat  is  supplied  to  the  node  by  direct  external  input  and  by  surface  reaction. 
External  input  may  take  several  forms,  such  as  frictional  heat  due  to  rubbing,  laser  irradiation, 
or  molten  metal  impingement.  External  heating  is  generally  considered  as  an  ignition  source; 
although  ignition  may  also  be  caused  by  aerodynamic  heating,  with  no  other  external  sources. 
Surface  reaction  refers  to  chemical  reactions  at  any  rate,  and  may  range  from  slow  oxidation  to 
very  rapid  combustion. 


For  the  second  case,  when  liquid  covers  the  surface,  Equation  (1)  is  written  in  two  distinct 
forms,  one  describing  the  energy  balance  at  the  liquid  surface,  the  other  describing  the  underlying 
solid  substrate.  In  the  liquid  layer  on  the  surface,  conduction  is  presumed  to  only  take  place 
between  the  liquid  and  the  solid  substrate.  Internal  energy  change  is  assumed  to  be  in  the  form 
of  phase  change,  and  the  energy  equation  takes  the  form: 


o ua 


dA  _ 

dt  qcond' 


surface  ’’’  Econv  + Eradiation  Einput  Ereaetion 


(2) 


The  thickness.  A,  is  considered  to  define  the  local  airfoil  thickness,  including  addition  or  losses 
due  to  recast  solidification  or  melting  of  the  surface,  and  H„(  is  the  latent  heat  of  fusion  of  the 
metal.  In  the  solid  substrate,  heat  flow  is  limited  to  thermal  conduction  between  the  nodal  mass, 
its  surrounding  nodes,  and  the  liquid  surface.  The  relevant  energy  equation  takes  the  form; 


dE 

Econd,  surface 


Econd,  nodes 


(3) 
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Equations  (1)  and  (3),  which  define  the  rate  of  sensible  heating  of  the  airfoil,  can  be  written 
in  terms  of  the  temperatures  of  the  nodal  masses,  thus  defining  the  airfoil  temperature  field. 
However,  the  resulting  set  of  equations  must  be  solved  simultaneously  with  Equation  (2),  which 
defines  the  local  rate  of  melting  or  solidification  at  the  surface,  since  heat  rejection  during  surface 
phase  change  provides  a major  source  of  energy  to  heat  the  substrate. 

The  various  heat  losses  and  gains  appearing  in  Equation  (2)  must  be  evaluated  in  terms  of 
the  local  environmental  and  airfoil  conditions.  Losses  by  convection  (qc)  and  radiation  (qR)  are 
predicted  by: 

qc  = hc  (T.-T,),  (4) 

and 


qr  = Fr<r(T,4  - Te4),  (5) 

where  hc  is  the  convective  film  coefficient,  T,  is  the  local  surface  temperature  and  Te  is  the 
surrounding  environmental  temperature,  which  is  approximated  by  the  gas  temperature,  T„.  The 
surface  temperature  in  the  region  of  liquid  metal  flow  is  assumed  equal  to  the  metal  melting 
temperature,  Tm.  The  quantity  Fr  is  the  overall  radiation  configuration  factor,  which  defines  the 
effect  of  geometry  and  properties  on  radiation.  The  local  temperature  field  in  the  solid  material 
is  transient  with  three-dimensional  heat  flow.  Assuming  spanwise  temperature  gradients  to  be 
small,  and  treating  the  temperature  field  as  quasi-steady,  the  differential  equation  for  thermal 
conduction  can  be  simplified  and  solved  to  yield  an  approximate  expression  for  conduction  from 
the  surface  into  the  solid  (qk).  The  resulting  heat  flow  is: 

q„  = -f-  (Tm  - Tn),  (6) 

where  Tn  is  the  nodal  temperature  in  the  underlying  solid,  combining  Equations  (2),  (4),  (5),  and 
(6): 


dA  _ h,.(Tm  - T,)  + Fr<T  (Tm4  - T,4)  - q,„pu,  - qrpBCtion 
dt  pHs, 


6k(Tm  - T„) 
pH„A 


As  molten  metal  flows  over  the  airfoil  surface,  heating  of  the  solid  material  by  conduction  from 
the  liquid  may  far  exceed  the  heating  rate  due  to  chordwise  conduction  within  the  airfoil.  The 
numerical  procedure  chosen  for  solution  of  the  airfoil  energy  equations,  Equations  (1)  and  (3), 
determines  the  temperature  field  at  discrete  times  during  the  process.  The  surface  energy 
equation.  Equation  (7),  is  solved  at  the  end  of  each  time  increment.  Since  the  surface  phase 
change  process  and  the  substrate  temperature  field  are  so  strongly  interdependent,  the 
alternating  (rather  than  simultaneous)  solution  of  the  energy  equations  results  in  a substantial 
error.  The  problem  is  alleviated  by  introducing  a subsidiary  energy  equation  for  the  nodal 
temperature  in  the  immediate  region  of  the  liquid  front  (downstream  edge  of  the  liquid  sheet), 
where  the  chordwise  temperature  gradient  is  most  severe,  and  where  the  transient  temperature 
rise  is  most  significant.  Equating  the  conduction  into  the  solid  to  the  rate  of  internal  energy 
change  in  the  node: 


T„). 


dT, 

dt 
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(8) 


Equations  (7)  and  (8)  are  solved  simultaneously  for  A and  Tn,  thus  defining  surface  phase  change 
and  the  change  in  nodal  temperature  in  the  immediate  vicinity  of  the  liquid  front;  however,  since 
these  equations  are  written  with  time  as  the  independent  variable,  an  auxiliary  equation  is 
required  to  define  the  motion  of  the  molten  metal  and  its  location  on  the  airfoil.  In  addition,  the 
convective  heat  transfer  and  chemical  reaction  at  the  surface  must  be  evaluated  prior  to  solving 
the  energy  equations. 


D.  AERODYNAMIC  PROCESSES 


Three  aerodynamic  processes  play  significant  roles  in  the  metal  combustion  phenomenon, 
namely  convective  heat  transfer  at  the  surface,  fluid  shear  stresses  and  drag  on  the  flowing  liquid 
metal,  and  diffusion  of  oxygen  to  the  metal  surface. 


Evaluation  of  these  processes  requires  an  evaluation  of  the  aerodynamic  boundary  layer  on 
the  surface.  The  approximate  methods  of  Walz  and  Buri  (Reference  6)  were  chosen  to  perform 
this  calculation,  yielding  satisfactorily  accurate  details  of  the  development  of  the  boundary  layer 
and  magnitude  of  the  shear  stress  acting  on  the  surface.  The  local  convective  heat  transfer  is 
determined  from  the  calculated  shear  stress  (friction  coefficient,  Cf)  by  use  of  the  Karman 
Boelter-Martinelli  analogy.  (Reference  7) 


PrRefl  V Cf/2 

5Pr  + 5 In  (5Pr  + 1)  + V2/C,  - 14 


(9) 


where  Nu0  = h0/k,  Re0  is  the  Reynolds  number,  Pr  is  the  Prandtl  number,  and  8 is  the  momentum 
thickness  evaluated  in  the  boundary  layer  analysis. 

Oxygen  diffusion  through  the  boundary  layer  is  determined  from:  (Reference  8) 

mc  = Kxm  ( * J (10) 

where  ri^  is  the  oxygen  mass  flux,  x„,.and  x,.*are  the  concentrations  at  the  wall  and  free  stream, 
and  Kxm  is  the  mass  transfer  coefficient  for  diffusion  determined  by  analogy  with  the  heat 
transfer  from: 


Kxm  = -gP°"gNu<9  ( -|1  )“  (11) 

Sc  and  Don  are  the  Schmidt  number  and  binary  diffusion  coefficient  for  oxygen-nitrogen 
mixtures. 


E.  REACTION  KINETICS 

The  reaction  rate  at  the  metal  surface,  either  as  slow  oxidation  or  as  rapid  combustion,  is 
governed  by  the  rate  at  which  oxygen  can  flow  from  the  air  to  the  surface  and  bond  with  free 
metal  in  the  substrate.  The  overall  process  is  governed  by  the  most  restrictive  step  in  the 
foregoing  chain  of  events.  The  most  restrictive  step  is  dependent  on  the  surface  temperature, 
amount  of  surface  oxides,  and  the  nature  of  the  boundary  layer  mass  diffusion  process. 
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The  governing  process  is  assumed  to  be  the  most  restrictive  of  the  following: 

1.  Oxygen  diffusion  through  the  surface  oxides  at  low  temperature  (parabolic 
kinetics),  dependent  on  surface  temperature  and  amount  of  accumulated 
oxides 

2.  Linear  kinetics  at  high  temperatures,  dependent  on  surface  temperature  only 

3.  Oxygen  diffusion  rate  through  the  aerodynamic  boundary  layer. 

The  kinetic  models  describing  these  three  conditions  are: 

1.  Parabolic  (Reference  9) 


29540 


dM  3.63  X 10’  e T 


dt  M 

2.  Linear  (Reference  10) 


, T<1300°K 


- 2.60  X 10>  e - , T>1300°K 


3.  Oxygen  Diffusion  in  Boundary  Layer 


= m„  (Eq  10) 


where  M = Oxygen  bonded  on  surface,  gm/cm1 
T = Surface  temperature,  °K 
t = Time,  sec. 

The  corresponding  reaction  rate  at  the  surface  is: 


(12) 


(13) 


qr  = H 


dM 

dt 


(14) 


where  H „„  i9  the  heat  of  reaction  per  mass  of  oxygen,  based  on  the  reaction  assumed  to  occur  at 
a particular  condition.  Post-test  analysis  of  burned  airfoil  specimens  indicated  the  primary 
reaction  product  resulting  from  slow  oxidation  is  TiO„  while  rapid  combustion  results  in 
approximately  equal  weights  of  TiO  and  Ti,0,.  The  corresponding  heats  of  reaction  have  been 
incorporated  in  the  model,  i.e.,  TiOa  is  assumed  to  result  from  solid  surface  oxidation,  while  50^ 
TiO-Ti,0,  is  assumed  to  result  when  combustion  occurs  on  a molten  surface. 


F.  LIQUID  METAL  FLOW 

Molten  metal,  driven  by  aerodynamic  forces,  is  observed  to  flow  over  the  surface 
downstream  of  the  airfoil  leading  edge.  A model  of  this  flow  was  developed,  based  on  the  liquid 
being  driven  only  by  aerodynamic  shear  stress  in  the  boundary  layer.  Based  on  the  liquid  metal 
viscosity  data  of  Magnitskiy,  (Reference  11)  several  cases  were  evaluated  and  compared  with  the 
experimental  data.  Predicted  velocities  of  the  liquid  front  were  found  to  be  low  by  several  orders 
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of  magnitude.  Re-examination  of  the  high  speed  movies  indicated  that  the  liquid  flow  is  rippled, 
and  fluid  motion  over  the  surface  is  in  the  form  of  unsteady  waves.  The  primary  driving  force  for 
this  type  of  motion  is  aerodynamic  profile  drag;  much  greater  than  the  corresponding  shear  stress. 

A new  model  is  formulated  in  which  liquid  moves  as  waves  over  the  wetted  surface,  as  shown 
in  Figure  19,  accumulating  at  the  downstream  edge  as  a liquid  front.  Motion  of  this  front  is 
controlled  by  a balance  of  aerodynamic  drag  due  to  the  airstream  and  surface  tension  forces  at  the 
interface  between  liquid  and  the  solid  surface.  As  the  front  moves  downstream  over  the 
nonwetted  surface,  solidification  of  part  of  the  liquid  film  occurs,  depending  on  the  relative  heat 
losses  by  conduction  and  convection,  and  the  heat  gain  from  reaction.  Motion  of  the  liquid  front 
is  intermittent,  liquid  accumulating  until  aerodynamic  forces  are  sufficient  to  exceed  surface 
forces,  at  which  time  the  liquid  front  moves  over  the  surface.  When  sufficient  mass  is  solidified 
to  reduce  the  size  of  the  liquid  wave  below  the  critical  size,  motion  of  the  front  ceases.  On  a time- 
averaged  basis,  the  motion  can  be  treated  as  steady,  being  controlled  by  the  rate  of  liquid  flow 
rate  reaching  the  front.  If  the  rate  of  liquid  melting  at  the  airfoil  leading  edge  is  m„  the  fraction 
lost  from  the  leading  edge  {,  and  the  local  airfoil  thickness  is  A,  then  the  local  rate  of  molten 
metal  flow  in  the  x-direction  along  the  surface  is: 

X 

m(x)  = (1  - ()  m,  - J p ^ dx  (15) 


Figure  19.  Liquid  Flow  Model 


Equation  (15)  is  an  expression  of  continuity  of  mass  along  the  surface.  At  the  edge  of  the  liquid 
front  x = xe,  and  the  liquid  metal  flow  rate  on  the  surface  is  m,. 

Assuming  that  the  liquid  front  is  configured  as  an  elongated  droplet,  as  shown  in  Figure  19, 
the  aerodynamic  force  necessary  to  overcome  surface  tension  force  and  produce  fluid  motion  is: 


25 


r 


where  A is  the  droplet  height,  f is  the  component  of  fluid  surface  tension, parallel  to  the  solid 
surface,  CD  is  the  aerodynamic  drag  coefficient,  and  p,  and  U„,  are  the  air  density  and  velocity. 
The  motion  of  the  liquid  front  is  related  to  its  thickness  and  flow  rate  by: 


m„  = p A V 


Solving  Eqs.  (15),  (16),  and  (17)  for  the  liquid  velocity: 


-{)“!-/  P dx 


Cd  Pa  U0 


which  defines  the  velocity  of  liquid  advance  over  the  surface. 


Q.  PROGRAM  LOGIC 


The  various  analytical  models  described  herein  have  been  incorporated  in  a computer 
program,  providing  an  integrated  procedure  for  the  prediction  of  ignition  and  self-sustained 
combustion  of  titanium.  The  operational  logic,  outlined  in  Figure  20,  allows  prediction  of  the 
ignition-combustion  sequence  for  a wide  range  of  ignition  sources  and  environmental  conditions. 


Set  Up  Airfoil  Geometry  and 
Grid  System 
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For  each  defined  set  of  environmental  and  ignition  conditions,  airfoil  geometry  and  material 
properties,  the  program  predicts  the  transient  temperature  field  during  the  duration  of  the 
ignition  source.  Ignition  is  presumed  to  have  occurred  if  the  transient  response  to  the  ignition 
source  results  in  a temperature  rise  to  the  melting  point  at  some  location  on  the  airfoil.  Self- 
sustained  combustion  is  presumed  if  reaction  heat  generation  is  sufficient  to  sustain  melting  after 
the  ignition  event  has  terminated.  If  ignition  and  self-sustained  combustion  is  found  to  occur,  the 
program  predicts  the  propagation  rate  as  combustion  spreads  across  the  airfoil,  until  the 
specimen  being  analyzed  is  totally  consumed,  or  extinction  of  burning  occurs. 

H.  ANALYTICAL  RESULTS 

Figure  21  shows  a comparison  of  the  basic  analytical  model  with  the  results  of  a combustion 
test  performed  in  a complementary  program  (Contract  N62269-76-C-0429).  The  test  specimen 
was  Ti-6A1-4V,  and  was  of  the  contoured  design  shown  in  Figure  6.  Although  chord  wise  burning 
rates  are  in  reasonable  agreement,  the  predicted  spanwise  propagation  rate  is  much  too  large.  The 
basic  analytical  model  predicts  the  rate  of  surface  melting  as  a function  of  the  relative  heat  gains 
and  losses  at  the  surface,  but  does  not  consider  the  direction  of  heat  flow  at  the  surface  and  within 
the  airfoil.  Hence,  the  basic  model  lacks  the  capability  to  predict  the  direction  of  leading  edge 
melting,  resulting  in  excessively  large  predicted  rates  of  spanwise  propagation  along  the  thin 
leading  edge.  This  discrepancy  indicates  the  need  for  refinement  of  the  model,  particularly  in  the 
areas  of  melt  flow  and  leading  edge  combustion. 
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Analytical  Simulation 


Test  Run  Results 


Time  = 0.15  sec 

t t 


Time  = 0.35  sec 


| = Direction  of  Air  Flow 
Specimen  Number  2-1 


Temperature 

Pressure 

Velocity 


= 429°C  (804°F) 
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Figure  21.  Comparison  of  Analytical  Simulation  With  Actual  Test  Results 
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8ECTI0N  V 


GENERALIZED  ANALYTICAL  MODEL 


To  achieve  the  needed  refinements  in  the  model,  the  analytical  study  was  extended  in 
duration  and  scope.  To  provide  a more  useful  predictive  model,  several  additional  goals  were 
added  to  the  program.  These  goals  were:  refinement  of  the  liquid  flow  model,  addition  of  alloy 
effects,  addition  of  accountability  for  aerodynamic  melt  loss  from  the  surface,  and  expansion  of 
the  model  to  account  for  centrifugal  body  forces  due  to  rotation.  These  goals  were  achieved  in  two 
steps;  development  of  the  model  refinements  for  a static  airfoil,  then  the  addition  to  account  for 
rotation.  However,  since  addition  of  rotational  effects  required  generalization  of  many  different 
areas  of  the  model,  only  the  final  resulting  model  will  be  described  in  this  report. 

A.  DESCRIPTION  OF  THE  MODEL 

Comparison  of  melt  flow  velocities,  evaluated  from  photographic  data,  with  thermal 
diffusion  rates  in  the  solid  airfoil  substrate  indicated  that  liquid  metal  flow  is  the  primary  mode 
of  thermal  propagation  in  a burning  airfoil.  For  purposes  of  comparison,  a simplified  one- 
dimensional thermal  analysis  was  conducted,  and  evaluated  over  the  range  of  observed  test  data. 
In  the  development  of  the  basic  analytical  model,  rippled  melt  flow  was  assumed  to  occur  over  the 
surface,  driven  by  aerodynamic  shear,  with  no  consideration  of  the  detailed  structure  of  such 
flow. 

In  the  generalized  model,  the  detailed  structure  of  the  melt  flow  is  examined,  and  a model 
developed  for  simultaneous  transport  of  momentum,  heat,  and  oxygen  species  from  the  liquid 
surface  to  the  liquid-solid  interface.  The  effects  of  oxides  on  the  relevant  physical  properties  are 
incorporated  in  the  model,  within  the  limitations  of  available  data.  Continuity  relations  are 
developed  in  integral  form  to  allow  evaluation  of  the  melt  flow  processes  over  the  surface  of  the 
airfoil,  Aerodynamic  melt  loss  considerations  are  included  in  the  continuity  equations,  and  a 
suggested  loss  model  is  proposed.  Finally,  the  model  geometry  is  expanded  to  incorporate  the 
effects  of  centrifugal  forces  generated  by  rotation,  as  in  the  case  of  compressor  rotors. 

All  description  of  the  model  development  will  be  based  on  the  inclusion  of  rotation  effect. 
Figure  22  illustrates  a typical  situation,  with  burning  initiated  at  the  leading  edge,  and  melt 
flowing  in  a compound  direction  toward  the  trailing  edge  tip.  This  figure  shows  the  salient 
features  of  the  process,  with  melt  streaming  from  the  leading  edge  of  the  burning  zone  to  the 
downstream  edge  of  the  melt  region,  henceforth  called  the  melt  leading  edge.  The  leading  edge  of 
the  burning  blade  will  be  called  the  stagnation  region. 

B.  THERMOPHYSICAL  PROPERTIES 

The  physical  properties  of  titanium  relevant  to  the  overall  combustion  model  include  the 
thermodynamic  properties,  defining  the  state  and  internal  energy,  and  the  transport  properties, 
defining  the  rate  of  the  various  transport  processes.  These  properties  in  general  are  dependent  on 
the  metal  temperature,  alloying  elements  and  amount  of  oxygen  present,  and  the  crystalline 
structure.  Because  of  insufficient  data,  the  present  analysis  has  not  included  alloying  effects,  but 
has  as  far  as  possible  included  the  effect  of  temperature,  phase,  and  oxygen  content.  The  effect 
of  temperature  and  oxygen  content  are  considered  in  evaluating  the  melt  viscosity  and 
thermodynamic  state,  and  the  effect  of  temperature  is  considered  on  the  thermal  conductivity 
and  oxygen  diffusivity. 
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Figure  23  shows  a phase  diagram  for  the  titanium-oxygen  system,  developed  by  Schofield 
and  Bacon  (Reference  12).  Oxygen  in  titanium  is  seen  to  act  er  ? phase  stabilizer  resulting  in  the 
occurrence  of  a-Ti  at  temperatures  above  the  normal  a-fi  transition,  and  increasing  the  liquidus 
temperature  significantly.  It  is  further  noted  that  the  presence  of  dissolved  oxygen  causes  the 
occurrence  of  a solidus-liquidus  temperature  range,  or  two-phase  region,  at  low  concentrations. 
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Figure  23  The  Ti-Q  Phase  Diagram  up  to  35  wt  rc  Oxygen  After  Schofield  and  Bacon 


A corresponding  thermodynamic  phase  diagram  was  constructed,  based  on  thermodynamic 
data  for  pure  titanium  (Reference  13). 


The  presence  of  oxygen  was  assumed  to  act  only  as  a phase-stabilizer  causing  the  effect  of 
the  energy  of  a-3  transition  to  be  shifted  to  the  melting  temperature  range.  While  the  resulting 
enthalpy  diagram,  shown  in  figure  24,  may  not  be  entirely  correct,  it  provides  a workable  basis  for 
modeling  the  heat  transfer  through  the  melt  flow. 


The  nature  of  the  two-phase  melting  region,  shown  in  Figure  23,  was  assumed  to  affect  the 
viscosity  of  the  melting  metal,  particularly  near  the  solid-liquid  interface.  The  data  of  Magnitskiy 
(Reference  11)  for  liquid  titanium  showed  the  viscosity  to  vary  inversely  with  temperature.  To 
extend  this  data  into  the  two-phase  region,  a hyperbolic  variation  of  viscosity  was  assumed  of  the 
form 


v (T,  <)>) 


I'm 


T - T, 


] 


, T,  < T < T 


(19) 


where  T,  and  T are  the  solidus  and  liquidus  temperatures  at  the  local  oxygen  concentration 
(4>),  and  vm  19  the  kinematic  viscosity  at  Tt.  For  temperatures  at  or  above  the  liquidus,  the  data 
of  Reference  11  is  considered  valid. 
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Data  was  not  available  on  the  effect  of  oxygen  on  the  other  transport  properties  in  the  liquid 
state,  nor  for  that  matter  for  pure  titanium  in  the  liquid  state.  The  data  of  Roe,  Palmer,  and  Opie 
(Reference  14)  for  oxygen  diffusivity  in  titanium  was  extrapolated  to  the  melting  temperature, 
and  treated  as  temperature  dependent  only.  The  data  of  Touloukian  (Reference  15)  for  thermal 
conductivity  was  treated  in  the  same  manner. 
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Figure  24.  Enthalpy  of  Titanium-Oxygen  System  as  Function  of  Temperature 

C.  STRUCTURE  OF  THE  MELT  FLOW 

In  the  flowing  melt  region  a liquid  layer  of  thickness  5 is  assumed  to  flow  on  the  surface, 
driven  by  aerodynamic  forces  in  the  direction  of  air  flow,  and  by  centrifugal  forces  in  the  radial 
outward  direction.  Oxygen  flows  in  the  aerodynamic  boundary  layer  to  the  liquid  surface,  reacts 
with  the  melt,  and  the  resulting  heat  and  oxides  diffuse  into  the  melt  layer.  Shear  stress  across 
the  melt  is  assumed  constant,  and  momentum  transport  is  considered  small  by  comparison. 
Figure  25  illustrates  the  features  of  the  flowing  liquid  layer.  The  coordinates  (x,  y,  z)  are  in  the 
chordwise,  normal,  and  spanwise  directions,  respectively.  The  velocity  components  (u,  v) 
represent  the  local  liquid  velocity  in  the  two  surface  directions  and  are  seen  to  vary  with  distance 
from  the  surface.  The  temperature  and  oxygen  concentration  at  the  outer  surface  are  Ti  and  0,, 
and  at  the  inner  surface  are  T0  and  0O.  If  the  outer  temperature  exceeds  the  liquidus  (superheated 
liquid),  the  saturated  liquid  state  is  assumed  to  occur  at  a distance  6i,  from  the  solid  interface, 
with  corresponding  temperature  and  concentration.  Oxygen  flows  into  the  melt  layer  at  a rate 
m„,  and  the  liquid-solid  interface  recedes  or  advances  (melting  or  solidification)  at  a velocity  U. 
If  melting  occurs  at  a single  temperature,  as  in  the  case  of  pure  oxygen-free  titanium,  it  is 
assumed  to  occur  at  the  liquid-solid  interface  with  a corresponding  discontinuity  in  the  local 
enthalpy  (latent  heat  of  fusion).  However,  if  a two-phase  melting  zone  exists  (solidus  and  liquidus 
not  equal),  the  enthalpy  is  assumed  to  increase  linearly  from  the  solid  to  the  liquid  state  over  the 
depth  of  the  melting  zone. 

The  nature  of  the  melt  surface  is  still  assumed  to  be  ripply,  the  amplitude  and  wavelength 
of  the  surface  waves  dependent  on  the  aerodynamic  shear  and  the  melt  thickness.  For  very  high 
shear  loading,  waves  will  crest  and  he  stripped  from  the  surface,  creating  spark  showers  as 
observed  in  the  tests. 
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Figure  25.  Local  Velocity  Profiles  and  Nomenclature  Used  in  Melt  Flow  Analysis 

D.  INTEGRAL  TRANSPORT  EQUATIONS 

In  Figure  22  the  melt  region  is  illustrated  with  a single  streamline.  The  melt  flow  field  can 
be  considered  as  composed  of  a set  of  streamlines  originating  in  the  stagnation  region  and  ending 
at  the  downstream  edge  of  the  melt  region,  as  shown  in  Figure  26.  If  a coordinate  system  is 
described  by  (S,  J),  where  S is  the  distance  along  the  streamline  and  J is  the  streamline  number, 
analysis  of  the  melt  flow  can  be  performed  along  each  streamline.  Let  the  liquid  flow  rate  in  a 
streamline  at  a point  (S,  J)  be  pP,  then  P is  a vector  quantity  aligned  in  the  direction  of  the 
streamline.  The  flow  vector  P can  be  described  as  composed  of  components  (f)  in  the  chordwise 
and  (h)  in  the  spanwise  directions.  Thus, 

P = Ff  + jh  (20) 

where  ! and  j are  unit  vectors  in  the  component  directions.  The  direction  of  the  streamline  can  be 
described  by  the  angle  (\p)  with  the  chord,  where 

tan  4,  = f/h  (21) 


The  melting  airfoil  adds  to  the  streamline  flow  in  the  amount  pU  per  unit  length  of 
streamline.  If  the  fraction  of  flow  rate  lost  to  the  airstream  is  £ per  unit  length,  then  the  rate  of 
change  of  flow  in  the  streamline  is 


r 


where  s is  a unit  vector  in  the  direction  of  the  streamline,  and  the  derivative  is  taken  in  the 
direction  of  flow.  This  expression  can  also  be  written  in  terms  of  the  flow  components  as: 


f -uw-ff 


^ = U sin  0 - £h 


(23a) 

(23b) 


Figure  26.  Streamline  Indexing  System 

In  a similiar  fashion,  the  flow  rate  of  oxygen  carried  along  the  surface  by  the  melt  is  (g).  and 
is  related  to  the  melt  flow  by 


g = pF<j> 


(24) 


where 


F = \fP  + h2  , 

and  0 is  the  local  bulk  concentration  of  oxygen  in  the  melt. 
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The  rate  of  change  of  g along  the  surface  is  given  by: 

"k'  = - £0,  pF  (25) 

The  concentration  at  the  liquid  surface  (0,)  is  presumed  to  describe  the  average  concentration  of 
the  melt  loss  at  the  surface. 

The  set  of  equations  (23a,  b)  and  (25),  describe  the  change  in  liquid  mass  flow  and  oxygen 
content  of  the  streamlines  as  they  flow  over  the  surface. 

E.  LOCAL  TRANSPORT  PROCESSES 

Integration  of  the  integral  transport  equations  requires  evaluation  of  the  initial  conditions, 
i.e.,  values  of  (f,  h,  g)  at  the  stagnation  point.  In  addition,  these  transport  functions  must  be 
related  to  the  local  transport  processes  to  allow  evaluation  of  the  physical  variables,  such  as  local 
velocity,  temperature,  and  concentration. 

1.  Axial  (Chordwiae)  Velocity  Profile 

The  melt  flow  was  assumed  to  be  laminar,  with  constant  shear  stress,  and  with  no  slip  at  the 
liquid-solid  interface.  Assuming  the  velocity  profile  to  be  represented  by  a second  degree 
equation: 

u = a + b»7  + ct7*  (26) 


where 


V = y/5. 

The  shear  stress  r,  must  be  equal  to  the  aerodynamic  shear,  determined  from  boundary 
layer  analysis,  and  is  related  to  the  velocity  profile  by: 

pv i du  I ...  pe„  du  I ^7) 

6 dy  y.j  6 dy  y.„ 


where  subscripts  T'  and  ‘o’  refer  to  evaluation  at  conditions  at  the  liquid  and  solid  surfaces.  With 
the  no-slip  conditions  (u  = o,  y = o)  and  Equation  (27),  the  constants  in  Equation  (26)  can  be 
determined,  yielding  the  velocity  equation: 

u(r,)  = (28) 

The  axial  component  of  flow  in  the  streamline  (per  unit  width)  is: 


6 

pf  = J pudv  = 6 

o 


r ptW 

■o  2m. 


dr? 


(29) 


Hence, 


f = 


T&’ 

6p, 


(30) 
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2.  Radial  Velocity  Component 


■ 


Fluid  motion  in  the  radial  direction  is  driven  by  centrifugal  body  force  (0),  and  resisted  by 
aerodynamic  shear  stress  (t0)  due  to  the  motion  of  the  liquid  relative  to  the  air  stream.  A simple 
analysis,  balancing  radial  forces,  shows  that  the  shear  stress  generated  by  body  forces  varies 
linearly  from  the  liquid  surface  to  a maximum  at  the  solid  surface.  The  aerodynamic  shear  stress 
is  found  to  be  constant  across  the  melt.  The  total  stress  is: 

t = 1 85  (1-t;)  - r0  (31) 


where 


0 — pRo>*  and  t0  = Cr  -^7^ 


(32) 


The  centrifugal  force  depends  on  the  radius  to  the  blade  mid-span  (R)  and  the  rotational 
speed  (u>)  of  the  compressor  stage.  The  aerodynamic  shear  stress,  is  a function  of  the  local  friction 
coefficient  (Cf)  and  the  relative  liquid  velocity  at  the  air-liquid  interface  (v,).  Approximating  the 
variation  of  viscosity  in  the  melt  by 


M = (33) 

the  laminar  relation  between  shear  stress  and  velocity  gradient  is  combined  with  Equation  (31) 
and  integrated  to  yield  the  velocity  profile: 


the  surface  velocity  (v,)  is  determined  by  combining  Equations  (32)  and  (34),  for  i\  = 1,  yielding: 


The  radial  flow  component  is  now  calculated  in  the  same  manner  as  the  axial  component, 
resulting  in  the  relation: 


h = 


M. 

12p, 


To  5* 

6pi 


(36) 


3.  Oxygen  Concentration  Profile 

The  distribution  of  oxygen  mass  fraction  (0)  is  assumed  of  the  same  form  as  the  velocity 
profile,  i.e.,  second  order  power  function  of  the  distance  from  the  solid  surface: 

d>(  r;)  = a + b)j  + ci7a 


36 


The  coefficients  are  evaluated  by  satisfying  the  following  conditions  at  the  solid  and  liquid 
surfaces: 

a.  Flow  rate  of  oxygen  into  the  liquid  surface  is  equal  to  the  diffusion  rate  into 
the  melt. 

b.  Diffusion  rates  in  the  liquid  and  solid  phases  are  equal  at  the  solid-liquid 
interface. 

c.  The  concentrations  in  the  liquid  and  solid  are  equal  at  the  interface. 

The  oxygen  diffusion  rate  into  the  solid  is  evaluated  by  solving  the  transient  diffusion 
(Laplace)  equations  for  the  case  of  a moving  boundary.  Two  dimensionless  groups  appear  in  the 
solution,  the  mass  transfer  Peclet  number  (A),  relating  the  mass  transport  by  convection  and 
diffusion,  and  the  transport  ratio  (D,  a form  of  mass  transfer  Reynolds  number.  These  quantities 
are  defined  by: 


A 


6 U 

D„ 


and 


(37) 


where  D„  and  D,  are  the  oxygen  diffusivities  in  titanium  at  the  temperatures  T„  and  T,. 

The  resulting  concentration  profile  must  be  related  to  a reference  value.  Relating  to  the 
liquid  surface: 


The  bulk  concentration,  (0),  is  obtained  by  integration  of  the  product  of  the  velocity  and 
concentration  profiles: 


/ 0^d» 1 

4>  = (39) 

/ Vdt, 


where  V is  the  liquid  velocity  vector.  Combining  Equations  (28),  (34),  and  (38),  and  performing 
the  integration  indicated  in  (39)  yields  a relationship  between  the  bulk  and  liquid  surface 
concentrations  (0,  0,).  The  surface  concentration  is  the  positive  root  of  the  following  quadratic 
equation: 


F,0i*  + l(2+A)  F,  - f F,]  0,  + 


( 1 + -f)'  (F.-01  + F,  - ( 


(40) 
= 0 
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Where: 


_ p5lW  v,! 

B=  r 'II=TT 


The  concentration  at  the  solid  interface  (ij=0)  is  readily  determined  from  Equation  (38)  as: 


<t>  o 


4.  Local  Heat  Transfer 


(41) 


Heat  transfer  into  the  flowing  melt  consists  of  the  sum  of  several  components,  heat 
generation  due  to  chemical  reaction  at  the  surface,  discrete  surface  heat  input  (i.e.,  laser 
radiation),  and  convective  and  radiant  exchange  with  the  surrounding  environment.  The 
components  for  reaction,  convection,  and  radiation  are  described  by  Equations  (14),  (4),  and  (5). 
Figure  27  shows  a typical  line  of  state  points  (dashed  curve)  across  a layer  of  melt.  Points  “0”  and 
“1”  correspond  to  the  solid  interface  and  outer  surface,  and  point  “i”  is  the  location  of  the 
saturated  liquid  state.  The  corresponding  distribution  of  thermodynamic  states  is  shown  in 
Figure  28.  The  physical  correspondence  of  these  state  points  is  shown  in  Figure  25.  In  the  region 
0 < y < j|,  melting  is  occurring  and  enthalpy  is  assumed  to  vary  linearly  with  distance  from  the 
solid-liquid  interface.  Thus, 


H - Ho  y 

H,  - Ho  5, 


(42) 


In  the  region  between  points  “i”  and  “1,”  melt  is  superheated  liquid.  Heat  flow  is  assumed  to  be 
constant  in  this  region,  thus  the  net  heat  input  at  the  outer  surface  is  equal  to  the  heat  flow  across 
the  surface  “i.”  The  heat  flow  in  the  melting  region  is  related  to  the  temperature  field  by 


k 


cPT 

dy* 


(43) 


which  describes  one-dimensional  conduction  with  a moving  boundary.  This  assumes  that 
conduction  normal  to  the  solid  surface  is  the  only  significant  heat  flow  in  the  melt.  A relationship 


38 


was  described  in  Equation  (7)  between  melting  rate  and  net  heat  flow  rate  in  the  melt.  This 
relationship  can  be  written  in  a generalized  vector  form  as: 


p H„  0 = q.  - qk, 


(44) 


where: 


q,  = net  surface  heat  flux  normal  to  surface 
q,,  = conduction  heat  flux  into  solid  substrate. 


Figure  27.  Typical  Temperature  — Concentration  Relationship  in  Melt 


For  melting  with  an  imposed  concentration  profile,  the  latent  heat  of  fusion 
becomes 


H„  = Hi  - Ho 


(45) 


Integration  of  Equation  (43),  in  combination  with  Equations  (42),  (44),  and  (45)  yields  a 
relationship  between  the  temperatures  T,  and  T0.  For  the  case  with  superheat: 


Qk 

2k 


(46) 


For  melting  without  superheat.  Equation  (46)  still  applies  with  5,  = &.  The  solid  interfacial 
temperature  T„  correspond  to  the  solidus  temperature  at  the  oxygen  content  of  the  interface. 
Because  the  surface  heat  flux  (q.)  is  a function  of  the  surface  temperature,  evaluation  of  Equation 
(46)  is  an  iterative  procedure.  Details  of  the  solution  are  shown  in  the  User’s  Manual. 
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Figure  28.  Typical  Temperature  — Enthalpy  Relationship  in  Melt 

F.  STAGNATION  REGION  ANALYSIS 


Integration  of  Equations  (23)  and  (25),  with  evaluations  of  the  corresponding  local 
processes,  defines  the  development  of  melt  flow  over  the  airfoil  surface.  However,  this  evaluation 
requires  that  conditions  at  the  stagnation  point  he  known  to  provide  initial  conditions  for  the 
integration.  Specifically,  three  quantities  are  required  to  start  the  calculation,  the  melt  thickness 
(5),  the  bulk  concentration  ($),  and  the  melting  velocity  at  the  leading  edge  (11).  Figure  29  shows 
a cross-sectional  view  of  the  leading  edge  with  the  relevant  features. 
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Figure  29.  Stagnation  Point  Region  of  Rurning  Airfoil 


Based  on  post-test  examination  of  burned  specimens,  the  cross-sectional  shape  of  a melting 
leading  edge  will  be  considered  an  elliptic  cylinder  with  an  aspect  ratio  (R  = width/length). 

The  initial  concentration  is  determined  by  expanding  Equation  (25),  taking  the  limit  as 
s— »o,  and  noting  that  the  aerodynamic  shear  stress  vanishes  at  the  stagnation  point.  This  results 
in  an  expression  for  the  bulk  stagnation  concentration: 


Consider  the  starting  point  for  streamline  (J  = l)  in  Figure  26,  noting  that  the  leading  edge 
is  yawed  to  the  direction  of  airflow.  Figure  30  shows  a larger  planform  view  of  the  leading  edge  and 
the  relative  directions  of  the  various  heat  fluxes  and  the  melting  velocity. 
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y = Angle  of  Attack 
0 = Leading  Edge  Melting  Velocity 


Figure  30.  Heat  Flaw  and  Velocity  Relationship  at  Burning  Leading  Edge 

The  leading  edge  is  at  a yaw  angle  (7)  to  the  airstream.  The  surface  heat  flux  (qs)  is  normal 
to  the  leading  edge,  and  the  internal  conduction  flux  (qk)  points  in  the  direction  of  maximum 
temperature  gradient  in  the  solid.  The  magnitude  and  direction  of  the  melting  velocity  is 
determined  by  the  vector  addition  defined  in  Equation  (44).  It  is  worthy  of  note  that  the  melting 
velocity  at  the  leading  edge  lies  in  the  plane  of  the  airfoil  leading  edge  at  an  angle  defined  by  the 
heat  transfer.  Over  the  remainder  of  the  airfoil  surface  the  melting  velocity  vector  is  oriented 
normal  to  the  airfoil  surface,  since  the  surface  and  internal  fluxes  are  aligned  everywhere  except 
at  the  leading  edge.  It  is  the  effect  of  the  internal  conduction  term  in  Equation  (44)  that  controls 
the  relative  burning  rates  in  the  chordwise  and  spanwise  directions,  observed  in  the  combustion 
tests. 


The  initial  melt  thickness  is  found  by  taking  Equations  (23a,  b)  to  the  limit  at  the 
stagnation  point.  Writing  these  equations  in  terms  of  components  normal  to  the  leading  edge,  and 
considering  the  variation  of  aerodynamic  shear  stress  and  centrifugal  force  as  the  stagnation  point 
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is  approached,  an  expression  is  derived  relating  initial  melt  thickness  to  the  other  parameters  at 
the  leading  edge.  The  resulting  equation  in  (6)  is: 


6"  1 + 62  Pa  U2  oo  dd^n  cos2  7 = 12  Ml  U (48) 

where  f)  is  defined  by  Equation  (32),  R and  A are  the  leading  edge  aspect  ratio  and  thickness.  The 
derivative  of  the  friction  coefficient  is  evaluated  at  the  stagnation  point  and  is  found  to  be  a 
function  of  R,  A,  and  the  leading  edge  Reynolds  number. 


G.  STREAMLINE  CONTINUITY 

Independent  analysis  of  flow  in  adjacent  streamlines  does  not  implicitly  avoid  continuity 
conflicts,  such  as  intersection  of  streamlines.  To  ensure  that  continuity  is  satisfied,  a condition  of 
streamwise  continuity  is  imposed  on  the  solutions  for  streamline  thickness.  If  the  width  of  a 
streamline  is  (W),  then  continuity  requires  that  (FW)  be  constant  along  the  streamline,  where  F 
is  the  magnitude  of  the  flow  vector.  This  statement  is  expressed  by: 


d 

ds 


W2 


P + h2  ) =0 


(49) 


If  local  rates  of  change  of  shear  stress  and  body  force  are  considered  secondary,  Equation  (49)  can 
be  expanded  into: 


do 

ds 


5(n  + h2) 


dW 

ds 


W(2P  -f  3h2  + hf  -^  ) 
r 


(50) 


This  relation  expresses  the  required  correction  in  the  calculated  growth  of  the  melt  thickness 
necessary  to  satisfy  streamline  continuity. 

H.  MELT  LEADING  EDGE  FLOW 


Figure  19  illustrated  the  main  features  of  the  liquid  flow  model  developed  as  part  of  the 
basic  analytical  model.  Modeling  of  the  melt  leading  edge  flow  was  based  on  a balance  of 
aerodynamic  drag  forces  and  surface  tension  between  the  melt  and  the  solid  wall. 

In  generalizing  the  model,  the  flow  becomes  two-dimensional  and  the  model  must  consider 
both  aerodynamic  and  body  forces.  Figure  31  shows  the  melt  leading  edge  region  isolated  as  a 
free-body.  The  sum  of  forces  acting  on  the  fluid  must  be  zero  for  steady  motion.  These  forces 
include  the  aerodynamic  drag  force  (FD),  the  component  of  surface  tension  parallel  to  the  solid 
surface  (f),  and  the  internal  forces  due  to  momentum  change  in  the  liquid  as  its  velocity  is 
reduced  within  the  leading  edge  region  (Flnt).  If  the  liquid  velocity  is  V,  as  it  approaches  the 
leading  edge  zone,  and  the  edge  moves  at  an  average  velocity  V,.,  the  internal  force  is  equal  to  the 
liquid  momentum  change,  or: 

F,nt  = m (V,  - V.)  = poV,  (V,  - V„).  (51) 

Considering  continuity  of  (low 

p^V,  = p AV,  = p F (52) 
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where  F is  the  magnitude  of  the  flow  vector  defined  in  Equation  (20).  Combining  these 
expressions, 

F.„,  = j£f-(y-1)  (53) 

The  aerodynamic  drag  is 

Fd  = CD  A cos2  0 (54) 


Summation  of  forces  yields  the  relation 


Figure  .V.  Leading  Edge  of  Melt 

This  is  solved  for  the  leading  edge  thickness.  A,  which  is  then  combined  with  Equation  (52) 
to  yield  V,,  the  velocity  at  which  melted  titanium  spreads  over  the  airfoil  in  the  direction  of  the 
streamlines.  In  Equation  (55),  the  flow  vector  is  evaluated  at  the  end  of  the  streamline.  Because 
the  end  of  the  streamline  is  advancing  with  time,  F is  changing  slightly  during  the  evaluation  of 
the  melt  motion.  The  analysis  nas  been  refined  to  account  for  this  effect,  and  the  procedure  is 
outlined  in  detail  in  the  User’s  Manual . 

I.  STAGNATION  POINT  DATA 

The  stagnation  region  analysis  of  the  generalized  model  has  been  used  to  evaluate  selected 
leading  edge  melting  rate  data  measured  in  the  test  program.  In  the  test  data  the  environmental 
conditions,  the  initial  geometry,  and  the  leading  edge  melting  velocity  were  known  or  measured 
directly.  The  instantaneous  leading  edge  geometry  was  estimated  from  the  photographic  data 
and  specimen  geometry.  The  leading  edge  aspect  ratio  (R)  was  assumed  to  he  0.5  for  this  analysis, 
which  appeared  to  result  in  reasonable  data  consistency.  The  analysis  yielded  the  melt  thickness 
at  the  leading  edge,  the  melt  surface  oxygen  concentration  (0,),  and  the  oxygen  diffusion  rate  into 
the  melt.  For  comparison,  the  oxygen  diffusion  model  for  the  aerodynamic  boundary  layer, 
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Equations  (10)  and  (11),  where  used  to  predict  the  oxygen  concentration  (mol-fraction)  in  the  air 
at  the  melt  surface. 

Figure  32  shows  the  predicted  melt  thickness,  dimensionalized  with  leading  edge  thickness 
and  correlated  with  leading  edge  aerodynamic  Reynolds  number.  Reynolds  number  is  a measure 
of  surface  shear  stress,  which  governs  the  thickness  of  melt  accumulation.  As  airstream  Reynolds 
number  increases,  the  “scrubbing”  of  the  surface  reduces  melt  thickness,  apparently  reaching  an 
asymptotic  limit. 


Reynolds  Number,  Thousands  K" 

Figure  .12  Relative  Thickness  of  Liquid  Metal  Film  at  Stagnation  Point 

Figure  33  compares  the  oxygen  concentration  at  the  air-liquid  metal  interface  as  a function 
of  Peclet  number  (A),  a measure  of  the  relative  magnitudes  of  melting  rate  and  diffusion  rate. 
The  oxygen  concentration  in  the  metal  is  seen  to  be  relatively  independent  of  melting  rate, 
whereas  the  concentration  in  the  air  boundary  layer  is  strongly  affected,  being  strongly  depressed 
with  increased  flow  rate  into  the  melt.  The  data  indicates  that  the  surface  oxidation  rate  is 
controlled  bv  diffusion  in  the  melt  at  low  melting  rates  (A  < 15),  and  by  diffusion  in  the 
aerodynamic  boundary  layer  at  higher  melting  rates.  This  requires  that  the  analytical  model 
consider  both  mechanisms,  and  select  the  controlling  model  for  evaluation  of  oxygen  How  into  the 
reaction. 


i 


Figure  34  shows  a correlation  of  the  predicted  oxygen  flow  in  the  melt  as  a function  of  the 
system  dynamics.  Oxygen  flow  is  correlated  by  the  Stanton  number  (St),  defined  by: 


St  = 


pU  ( <t> , - 0) 


(56) 


where  oxygen  flow,  melting  velocity,  and  the  oxygen  concentrations  are  all  evaluated  at  the 
stagnation  point.  The  data  shown  in  the  figure  is  well  represented  by  the  empirical  correlation: 


St  = 7.17  A 


+ 0.00673A 


(57) 


! 
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Although  the  curve  is  predicted,  the  data  points  are  shown  to  indicate  the  range  of  test  data.  The 
relationship  shown  in  the  figure  between  Stanton  number  and  Peclet  number  takes  the  same  form 
as  the  data  in  the  published  literature  for  heat  transfer.  Heat  transfer  data  in  the  literature,  for 
stagnation  point  flow  and  for  boundary  layer  flow,  were  examined  and  found  to  be  correlated  by 
a single  empirical  relation  of  the  same  form  as  Equation  (57).  Based  on  this  result,  extension  of 
the  range  of  applicability  of  Equation  (57)  beyond  the  stagnation  point  is  _,nsidered  justified. 

J.  AIRFOIL  HEAT  TRANSFER 

The  energy  equation,  Equation  (1),  for  conduction  of  heat  within  the  solid  airfoil  has  been 
re-written  in  finite  difference  form  to  include  the  heat  input  from  the  flowing  melt,  defined  by 
Equation  (6).  A stability  analysis  to  determine  criteria  for  computational  stability  indicated  the 
most  stringent  condition  to  be  imposed  when  the  surface  is  covered  by  molten  metal.  Details  of 
the  finite  difference  equations  and  stability  criteria  are  shown  in  the  User’s  Manual. 

K.  COMPUTING  PROCEDURES 

Computing  procedures  in  the  generalized  model  have  been  modified  somewhat  from  those 
outlined  in  the  basic  model.  Procedural  details  are  shown  in  the  User’s  Manual,  but  an  outline 
of  the  logic  flow  is  shown  in  Figure  35.  In  general,  the  program  proceeds  through  several  main 
functional  areas  for  each  time  increment  in  the  simulation. 

These  function  areas  are  as  follows: 

1.  Evaluation  of  temperature  field  in  solid  airfoil, 

2.  test  for  occurrence  of  ignition, 

3.  analysis  of  flowing  melt  processes, 

4 evaluation  of  spread  of  melt-covered  region, 

5.  evaluation  of  change  in  airfoil  geometry  due  to  melting,  and 

6.  return  to  (I). 

Figure  35  shows  the  operational  flow  and  required  conditional  transfers.  Logic  diagrams  for 
the  various  operational  blocks  are  expanded  in  Ref.  16. 
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SECTION  VI 


SUMMARY  OF  RESULTS 


The  self-sustained  combustion  of  a selected  titanium  alloy  (Ti-8Al-lMo-lV),  has  been 
investigated  over  a range  of  temperatures,  pressures  and  flow-velocities  encompassing  the  typical 
gas  turbine  engine  compressor  environment.  Based  on  this  experimental  and  analytical  effort,  the 
interrelationship  of  the  several  important  environmental  parameters  may  be  summarized  as 
follows: 

1.  The  total  temperature  and  the  Reynolds  number  at  the  airfoil  leading  edge 
appear  to  be  satisfactory  correlation  parameters  for  defining  the  lower  limit 
of  self-sustained  combustion. 

2.  Burning  rate  in  the  chordwise  direction  is  strongly  dependant  on  Reynolds 
number. 

3.  Total  damage  is  a function  of  spanwise  burning  rate,  which  in  turn  is  a 
complex  function  of  pressure,  velocity,  and  temperature.  Damage  does  not 
appear  to  be  correlated  by  Reynolds  number,  but  generally  increases  with 
temperature,  an  effect  of  the  reduced  heat  sink  cooling  capacity  of  the  airfoil. 

4.  Examination  of  high-speed  film  indicates  the  importance  of  melt  removal  by 
aerodynamic  forces.  At  low  velocities  and  reduced  combustion  rates, 
considerably  quantities  of  liquid  are  retained  and  flow  over  the  specimen 
promoting  continued  combustion,  and  are  removed  in  relatively  large 
globules.  With  increasing  velocity,  greater  aerodynamic  forces  promote  more 
rapid  removal  of  burning  droplets,  which  appear  to  be  of  progressively  finer 
size.  This  behavior  has  considerable  impact  on  the  extent  of  downstream 
damage  (fresh  ignition  sites)  in  a closed  environment  such  as  a gas  turbine 
compressor. 

5.  Comparative  tests  performed  at  the  same  environmental  conditions,  but 
with  different  orientation  of  the  specimens,  show  very  significant  effects  of 
the  aerodynamic  angle-of-attack  on  the  nature  of  the  propagation.  Specimen 
damage  appears  greater  with  angled  than  with  parallel  airflow.  This  effect 
appears  to  be  due  to  changes  in  the  nature  of  the  airflow  field  around  the 
specimen,  possibly  with  separation  or  secondary  flows  in  the  spanwise 
direction,  which  results  in  rapid  spanwise  spreading  of  liquid  metal  over  the 
surface.  Although  the  result  on  combustion  is  apparent,  the  cause-effect 
relation  is  not  fully  understood  at  this  time. 

A basic  analytical  model  has  been  developed  to  simulate  the  ignition  and  self-sustained 
combustion  of  stationary  titanium  alloy  airfoils  in  gas  turbine  environments,  and  the  model  has 
been  programmed  for  computer  evaluation.  Preliminary  results  with  the  model  indicate  that  it  is 
generally  in  good  agreement  with  the  test  data  in  predicting  the  lower  limits  of  self-sustained 
combustion  and  in  simulating  liquid  metal  flow;  however,  some  details  of  the  model  required 
further  development  to  provide  a predictive  technique  generally  applicable  to  a wide  range  of  gas 
turbine  engine  environments. 

The  analytical  model  was  further  generalized  and  programmed  for  computer  solution  to 
include  the  effects  of  centrifugal  body  forces  on  liquid  metal  flow  and  propagation,  thus  extending 
its  application  to  analysis  of  compressor  stators  and  rotors.  Also,  to  provide  better  agreement  with 


observed  experimental  behavior,  liquid  flow  models  were  expanded  to  include  aerodynamic 
coupling  mechanisms  controlling  loss  of  melt  from  the  surface,  and  physical  property  models  were 
developed  to  include  the  effect  of  absorbed  oxygen  on  thermodynamic  and  transport  properties. 
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Notation 

CD 

CF 

Cp 

D 

Don 

E 

P 

F 

Fr 

f 

g 

h 

hr 

H 

H,„o 

H„ 


i.j 


k 

K»m 

M 

m 

m„ 


m, 


n 


= aerodynamic  drag  coefficient 
= coefficient  of  friction 
= specific  heat 
= oxygen  diffusivity  in  metal 
= binary  oxygen-nitrogen  diffusion  coefficient 
= internal  energy 
- melt  flow  vector 
= melt  flow  scalar 
= radiation  configuration  factor 
= chordwise  component  of  melt  flow 
= oxygen  convective  flow  in  melt 
= radial  component  of  melt  flow 
= convective  heat  transfer  coefficient 
= enthalpy  of  metal 

= heat  of  reaction  per  unit  mass  of  oxygen 
= latent  heat  of  fusion 

unit  vectors  in  chordwise  and  radial  directions 
= thermal  conductivity 
= mass  transfer  coefficient 
= mass  of  oxygen  bonded  on  surface 
= mass  flow  rate 
= mass  flow  rate  at  liquid  edge 
= mass  flow  rate  at  blade  leading  edge 
= oxygen  mass  flux 
= direction  normal  to  surface 
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LI8T  OF  8YMBOL8 

Notation 


A»'  Q» 

= 

surface  heat  flux  vector,  scalar 

Qk.  Qk 

= 

internal  conduction  vector,  scalar 

R 

= 

leading  edge  aspect  ratio 

R 

= 

radius  to  mid-span  of  blade 

S,J 

= 

streamline  position  and  location  indices 

s 

= 

distance  along  streamline 

t 

= 

time 

T 

= 

temperature  in  melted  metal 

T 

1 k 

= 

gas  temperature 

T 

1 m 

= 

melting  temperature 

Tn 

= 

node  temperature 

U 

= 

melting  velocity  of  surface 

= 

free-stream  air  velocity 

U,  V 

= 

chordwise  and  radial  velocity  components  (melt) 

v,v 

= 

liquid  metal  velocity  vector,  scalar 

W 

= 

width  of  streamline 

X.Y.Z 

= 

chordwise,  normal,  radial  coordinate  system 

X, 

= 

end  of  liquid  metal  region 

n 

= 

centrifugal  body  force  per  volume 

y 

= 

leading  edge  yaw  angle 

V 

= 

mass  transport  ratio  (Eq.  37) 

<5 

= 

liquid  metal  film  thickness 

fh 

= 

thickness  of  two-phase  region 

A 

= 

local  thickness  of  airfoil 
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To 
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X<„  „ 


X„, 


on 
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u > 

u 

Mn 


N„ 


Pr 

Rew 


Sc 


St 


= component  of  surface  tension  parallel  to  surface 
= dimensionless  distance  from  solid  surface 
= boundary  layer  momentum  thickness 
= thickness  of  liquid  metal  leading  edge 
= mass  transfer  Peclet  number 
= dynamic  viscosity  of  liquid  metal 
= kinematic  viscosity  of  liquid  metal 
= kinematic  viscosity  at  liquidus  temperature 
= fraction  of  melt  flow  loss  from  surface 
= metal  density 
= air  density 

= Stefan-Boltzmann  constant 
= aerodynamic  shear  stress  in  chordwise  direction 
= aerodynamic  shear  stress  in  radial  direction 
= local  oxygen  mass-fraction  in  liquid  metal 
= bulk  (mean)  oxygen  mass-fraction  in  liquid  metal 
= oxygen  mol-fraction  in  air  at  surface 
= oxygen  mol-fraction  in  air  free-stream 
= streamline  angle  relative  to  chord 
= rotational  speed  of  blade 
= angle  of  attack  of  test  specimen 
= Mach  number 

= Nusselt  number  for  heat  transfer 

Prandtl  number 

Reynolds  number  based  on  momentum  thickness 

= Schmidt  number 

= Stanton  number  for  mass  transport 
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Notation 


Subscripts 


interface  of  liquid  and  two-phase  regions 

solid 

liquid 

liquid-solid  interface 
liquid  metal-air  interface 
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SECTION  I 


INTRODUCTION 


Prfltt  & Whitney  Aircraft  Group  Customer  Computer  Deck  (CCD1152-0.0)  is  an  analytical 
model  for  the  simulation  of  ignition  and  self-sustained  combustion  of  titanium  and  its  alloys  in 
gas  turbine  engine  compressor  environments. 

The  object  of  the  model  is  to  predict  the  susceptibility  of  titanium  to  be  ignited  and  undergo 
self-sustained  combustion.  In  addition,  if  combustion  occurs,  the  model  predicts  the  character- 
istics and  rate  of  propagation  across  the  burning  surfaces. 


The  program  provides  the  capability  of  analyzing  either  compressor  rotors  or  stators,  with 
considerable  flexibility  in  varying  the  geometry,  material,  and  environmental  conditions.  Several 
different  ignition  modes  can  be  considered,  including  external  heat  input,  aerodynamic  heating, 
mechanical  rubbing,  fracture,  or  melt  impingement.  Location  of  ignition  can  be  specified  or,  in 
the  case  of  aerodynamic  heating,  the  model  determines  the  location. 

The  User’s  Manual  describes  the  Titanium  Airfoil  Combustion  Program,  its  functional 
characteristics  and  operation,  and  the  diagnostics  necessary  to  assist  in  installation  and  checkout 
of  the  program  in  the  user’s  facility. 


SECTION  II 


DESCRIPTION  OF  THE  PHYSICAL  SYSTEM 


Titanium  is  observed  to  burn  in  self-sustained  combustion  when  ignited  under  certain 
environmental  conditions.  In  particular,  the  conditions  occurring  in  current  axial-flow  gas 
turbine  engine  compressors  have  been  found  to  be  favorable  for  the  propagation  of  titanium  fires. 
Figure  II-l  illustrates  the  features  of  a burning  compressor  blade,  ignited  at  the  midspan  leading 
edge.  In  this  schematic,  the  leading  edge  of  the  blade  is  melting,  and  the  liquid  metal  is  driven 
by  aerodynamic  and  centrifugal  forces  over  the  surface  of  the  blade.  Burning  takes  place  at  the 
liquid  surface,  generating  heat  which  causes  melting  of  the  airfoil.  The  primary  process  for  fire 
propagation  is  the  flow  of  burning  liquid  metal.  The  rate  and  direction  of  spreading  is  defined  by 
the  combined  thermal  and  mass  transfer  processes.  The  thermal  process  is  controlled  by  the 
balance  of  heat  gains  (heat  input  and  reaction)  and  heat  losses  (conduction,  convection,  and 
radiation),  which  together  define  the  rate  of  melting  Two  mass  transfer  processes  occur,  the 
transport  of  oxygen  from  the  air  into  the  liquid  metal,  controlling  the  reaction  rate,  and  the  flow 
of  liquid  metal  over  the  airfoil  surface,  controlling  the  rate  of  propagation. 

Melting  occurs  at  the  liquid-solid  interface,  and  the  resulting  liquid  spreads  across  the 
airfoil  along  the  “melt  edge,”  as  shown  in  Figure  II- 1 . In  the  case  of  nonrotating  stators,  liquid 
flow  is  in  the  chordwise  direction  only,  driven  by  aerodynamic  friction.  Liquid  metal,  stripped 
from  the  surface  by  the  airstream,  or  streaming  from  the  trailing  edge  of  the  burning  airfoil, 
provides  the  energy  source  to  cause  ignition  of  subsequent  downstream  compressor  stages. 


SECTION  III 


THE  ANALYTICAL  MODEL 


The  analytical  model  of  the  physical  system  consists  of  several  parts,  providing  a 
simultaneous  analysis  of  the  several  mechanisms  controlling  the  combustion  process.  These 
several  parts  describe  the  aerodynamic  processes  external  to  the  burning  airfoil,  the  reaction 
kinetics  of  the  burning  metal,  heat  transfer  within  the  airfoil,  melting  of  the  airfoil  and  flow  of  the 
molten  metal,  diffusion  of  oxygen  into  the  metal,  and  changes  in  the  airfoil  geometry  due  to 
melting. 

Models  are  written  in  finite  difference  form,  and  simulation  of  the  process  procedes  as  a 
marching  solution  in  time  and  space.  Iterative  procedures  are  used  in  several  subsystems 
requiring  simultaneous  solutions  of  implicit  equations. 

Development  of  the  analytical  models  are  described  in  Reference  1 and  the  descriptive 
equations  and  solution  logic  are  outlined  in  Appendices  A and  B. 
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SECTION  IV 


DESCRIPTION  OF  THE  COMPUTER  PROGRAM 


A.  PURPOSE 

The  purpose  of  the  computer  program  is  to  provide  an  analytical  technique  for  evaluating 
the  susceptibility  of  titanium  alloys  to  ignite  and  burn  when  used  in  gas  turbine  engine 
components.  Specifically,  the  program  is  designed  for  analysis  of  axial  flow  compressor  blades 
and  vanes  in  the  aerodynamic  environment  of  high  performance  GTE  compressors. 

B.  GENERAL  DESCRIPTION 

The  overall  system  logic  is  shown  in  the  flow  diagram  of  Figure  IV-1 . Detailed  flow  diagrams 
of  the  major  subsystems  are  shown  in  Appendix  A.  The  program  operates  in  two  major  sections, 
(a)  evaluation  of  the  temperature  distribution  in  the  solid  metal  and  (b),  evaluation  of  the 
melting  and  associated  liquid  metal  processes.  The  first  section  operates  in  a time  stepping  mode, 
whether  or  not  burning  is  occurring.  If  burning  occurs,  the  liquid  metal  section  operates  for  each 
time  step  to  evaluate  the  melting  process  and  spread  of  molten  metal,  and  provides  input  for 
evaluation  of  the  temperature  field  at  the  next  time  increment. 

C.  FUNCTIONS  PERFORMED 

Referring  to  Figure  IV-1,  the  sequence  of  functional  operation  is  in  the  following  order: 

1.  Input  geometry,  environmental  conditions,  material  properties,  ignition 
mode  and  energy  level, 

2.  Initialize  systems:  set  nodal  temperatures,  velocity  profile, 

3.  Calculate  geometric  quantities  (nodal  masses,  areas,  conduction  path 
lengths), 

4.  Calculate  aerodynamic  boundary  layer,  shear  stresses,  heat  and  mass 
transfer  coefficients, 

5.  Calculate  surface  heat  generation, 

6.  Evaluate  temperature  field  in  airfoil, 

7.  Test  for  melting  (if  none,  return  to  (3)), 

8.  If  melting  (and  burning)  occurs,  perform  leading  edge  (stagnation  point) 
analysis-calculate  liquid  metal  film  thickness,  oxygen  absorption  rates  and 
diffusion  profiles,  melting  rate,  temperature  profile, 

9.  Calculate  same  quantities  along  liquid  flow  streamlines, 

10.  Calculate  spread  of  burning  liquid  over  airfoil  surface, 

11.  Some  liquid  metal  removed  from  surface  by  air  stream, 

12.  Evaluate  change  in  airfoil  geometry  due  to  melting,  and 

13.  Increment  time  and  return  to  (3). 
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Figure  IV-1.  Program  Logic  Diagram 
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The  procedure  continues  until  a defined  total  simulation  time  is  reached.  In  the  flow 
diagram  (X,Z)  are  coordinates  in  the  chordwise  and  radial  (spanwise)  directions.  Indices  (J,  S) 
refer  to  streamline  index  and  incremental  distance  along  streamline  from  airfoil  leading  edge,  as 
shown  in  Figure  IV-2.  The  subscripted  indices  (Je,  Se),  (Xe,  Ze)  refer  to  quantities  at  the  “melt 
edge.” 

Detailed  description  of  the  functions  and  functional  operation  are  presented  in  Appendix  A. 


Figure  IV-2.  Streamline  Indexing  System 


SECTION  V 


GEOMETRIC  NOMENCLATURE 


Figure  V-l  illustrates  the  nodal  breakup  of  an  airfoil,  and  the  significant  geometric 
parameters.  The  number  of  spatial  increments  is  currently  fixed  in  the  program,  as  shown  in 
Figure  V-l. 
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Figure  V-l.  Nomenclature  of  Airfoil  Geometry 
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SECTION  VI 
INPUT/OUTPUT 

The  program  requires  as  input  the  specification  of  airfoil  geometry,  material  properties, 
environmental  conditions,  and  mode  of  ignition  energy  input.  The  following  input  description 
illustrates  the  required  format  for  card  input  of  the  specifying  parameters. 

Program  output  provides  graphic  and  tabulated  data  on  ignition  occurrence,  subsequent 
rates  of  airfoil  damage  and  fire  propagation,  and  resulting  radial  and  axial  carry-over  of 
combustion  products  throughout  the  engine.  Approximate  graphic  display  is  plotted  on  the 
printed  output,  and  accurate  graphic  data  is  optionally  available  through  the  use  of  a Calcomp 
plotter. 

A.  INPUT  DESCRIPTION 


Card 

Columns 

Variable 

Definition 

1 

1-80 

TITLE 

Any  descriptive  information  to  identify  this  case. 

2 

1-10 

STPTIM 

Total  simulation  time  in  seconds. 

11-20 

TINC 

Time  increment  for  graphic  output  in  seconds. 

3 

1-10 

DTL 

X-direction  distance  from  the  leading  edge  of  the  root 
to  the  leading  edge  of  the  tip  in  inches  (sweep  length). 
See  Figure  V-l. 

11-20 

AVSP 

Average  airfoil  span  in  inches. 

21-30 

TC 

Tip  chord  in  inches. 

31-40 

RC 

Root  chord  in  inches. 

41-50 

ARAT 

Aspect  ratio  of  blade. 

51-60 

RDISC 

Distance  from  center  of  rotation  to  blade  root  in 
inches. 

61-70 

VTAN 

Tangential  velocity  of  blade  in  ft/sec. 

71-80 

TSTRT 

Melting  temperature  of  the  metal  in  °R. 

4 

1-10 

TCOOL 

Free  stream  gas  temperature  in  °R. 

11-20 

PRES 

Free  stream  gas  pressure  in  psia. 

21-30 

THOXID 

Initial  oxide  thickness  in  inches. 

31-40 

RHOM 

Density  of  solid  titanium  in  (tbm/in!). 

41-50 

RHOL 

Density  of  liquid  titanium  in  (tbm/in!). 

51-60 

KUQ 

Thermal  conductivity  of  liquid  titanium  in  (BTU/hr- 
ft-°R). 

61-70 

TSINK 

Temperature  to  which  blade  is  radiating  in  °R.  If 
input  equal  to  zero,  the  gas  temperature  will  be  used. 

71-80 

FEFA 

Product  of  view  factor  and  emissivity  to  be  used  for 
radiation  heat  transfer  between  blade  and  gas  en- 
vironment (Fr). 

5 

1-10 

STEN 

Surface  tension  in  (tbf/ft). 

11-20 

FLOS 

Fraction  of  melt  flow  lost  from  surface  (%  loss/ft). 

21-30 

CD 

Coefficient  of  drag. 
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B.  METAL  THERMAL  CONDUCTIVITY  CURVE 


Card  Columns 

Variable 

Definition 

6 1-2 

IK 

Number  of  points  in  metal  conductivity  curve  (max- 
imum of  20  points,  right-adjusted). 

3-10 

VALKPT(l) 

First  temperature  point  in  curve  in  °R. 

11-20 

VALK(l) 

Conductivity  corresponding  to  first  temperature  in 
(Btu/hr-ft-°R). 

21-30 

VALKPT(2) 

Second  temperature  point  in  curve. 

31-40 

VALK(2) 

Conductivity  corresponding  to  second  temperature. 

41-50 

VALKPT(3) 

Same  as  VALKPT(2). 

51-60 

VALK(3) 

Same  as  VALK(2). 

61-70 

VALKPT(4) 

Same  as  VALKPT(2). 

71-80 

VALK(4) 

Same  as  VALK(2). 

If  IK  >4,  continue  on  another  card  in  fields  of  10.  Omit  input  of  IK  on  successive  cards. 

C.  METAL  SPECIFIC  HEAT  CURVE 

Card 

Columns 

Variable 

Definition 

7 

1-2 

ICP 

Number  of  points  in  metal  specific  heat  curve 
(maximum  of  20  points,  right -adjusted). 

3-10 

CPTVAL(l) 

First  temperature  point  in  curve  in  °R. 

11-20 

CFVAL(l) 

Specific  heat  corresponding  to  first  temperature  in 
(Btu/tbm-°R). 

21-30 

CPTVAL(2) 

Second  temperature  point  in  curve. 

31-40 

CPVAL(2) 

Specific  heat  corresponding  u.  second  temperature. 

41-50 

CPTVAL(3) 

Same  as  CPTVAL(2). 

51-60 

CPVAL(3) 

Same  as  CPVAL(2). 

61-70 

CPTVAL(4) 

Same  as  CPTVAL(2). 

71-80 

CPVAL(4) 

Same  as  CPVAL(2). 

If  ICP  >4,  continue  on  another  card  in  fields  of  10.  Omit  input  of  ICP  on  successive  cards. 


D.  TITANIUM  LIQUID  AND  SOLIDUS  TEMPERATURE  CURVE 


Card  Columns 

Variable 

Definition 

8 1-2 

NTS 

Number  of  points  in  solid  and  liquid  temperatures 
as  a function  of  oxygen  concentration  (maximum  of 
20  points,  right-adjusted). 

3-10 

PHIV(l) 

First  oxygen  mass  fraction  point  in  curve. 

11-20 

TLIQ(l) 

Temperature  of  liquid  corresponding  to  first  ox- 
ygen mass  fraction  in  °R. 

21-30 

TSOL(  1 ) 

Temperature  of  solid  corresponding  to  first  oxygen 
mass  fraction  in  °R. 

31-40 

PHIV(2) 

Second  oxygen  mass  fraction  point  in  curve. 

41-50 

TLIQI2) 

Temperature  of  liquid  corresponding  to  second 
point  in  curve. 

51-60 

TSOL(2) 

Temperature  of  solid  corresponding  to  second  point 
in  curve. 

If  NTS  >2,  continue  on  another  card  in 
sets  of  input  per  card. 

fields  of  10.  Omit  input  of  NTS  on  successive  cards.  Two 

E.  ROOT  THICKNESS  VARIATIONS 

Card 

Columns 

Variable 

Definition 

9 

1-10 

TCORD(l,l) 

Root  thickness  at  nodal  location  1 in  inches 

11-20 

TCORD(2,l) 

Root  thickness  at  nodal  location  2. 

21-30 

TCORD(3,l) 

31-40 

TCORD(4,l) 

41-50 

TCORD(5,l) 

51-60 

TCORD(6,l) 

61-70 

TCORD(7,l) 

71-80 

TCORD(8,l ) 

Root  thickness  at  nodal  location  8. 

10 

1-10 

TCORD(9,l) 

Root  thickness  at  nodal  location  9. 

11-20 

TCORD(10,1) 

21-30 

TCORD(ll.l) 

Root  thickness  at  nodal  location  11. 

F.  TIP  THICKNESS  VARIATIONS 


Card 

Columns 

Variable 

Definition 

11 

1-10 

TCORD(l,2) 

Tip  thickness  at  nodal  location  1 in  inches. 

11-20 

TCORD(2,2) 

Tip  thickness  at  nodal  location  2. 

21-30 

TCORD(3,2) 

31-40 

TCORD(4,2) 

41-50 

TCORD(5,2) 

51-60 

TCORD(6,2) 

61-70 

TCORD(7,2) 

71-80 

TCORD(8,2) 

12 

1-10 

TCORD(9,2) 

Tip  thickness  at  nodal  location  9. 

11-20 

TCORD(10,2) 

21-30 

TCORD(  11,2) 

Tip  thickness  at  nodal  location  11. 

G.  ROOT  VELOCITY  PROFILE 


Card 

Columns 

Variable 

Definition 

13 

1-10 

U(l,l) 

Root  velocity  at  nodal  location  1 in  (ft/sec) 

11-20 

U(2,l) 

Root  velocity  at  nodal  location  2. 

21-30 

U(3,l) 

31-40 

U(4,l) 

41-50 

U(5,l) 

51-60 

U(6,l) 

61-70 

U(7,l) 

71-80 

U(8,l) 

Root  velocity  at  nodal  location  8. 

14 

1-10 

U (9, 1 ) 

Root  velocity  at  nodal  location  9. 

11-20 

U(10,l) 

21-30 

U(ll,l) 

Root  velocity  at  nodal  location  11. 

H.  TIP  VELOCITY  PROFILE 


Card 

Columns 

Variable 

Definition 

15 

1-10 

U(l,2) 

Tip  velocity  at  nodal  location  1 in  (ft/sec). 

11-20 

U(2,2) 

Tip  velocity  at  nodal  location  2. 

21-30 

U(3,2) 

31-40 

U(4,2) 

41-50 

U(5,2) 

51-60 

U(6,2) 

61-70 

U(7,2) 

71-80 

U(8,2) 

16 

1-10 

U(9,2) 

Tip  velocity  at  nodal  location  9. 

11-20 

U(10,2) 

21-30 

U(1 1,2) 

Tip  velocity  at  nodal  location  11. 

16 

1-5 

IHEAT 

Type  of  heat  input. 

IHEAT  = 1,  laser  input. 

IHEAT=2,  aerodynamic  input. 

IHEAT=3,  rub. 

IHEAT=4,  metal  impingement. 

IHEAT =5,  fracture. 

Value  must  be  a right -adjusted  integer. 

6-10 

NODES 

Number  of  nodes  which  are  effected  by  heat  input. 

Not  required  if  IHEAT=2.  Value  must  be  a right- 
adjusted  integer. 

11-20 

DTGAS 

Input  if  IHEAT  = 2.  Increase  in  gas  temperature  °R. 

21-30 

QHIN 

Input  if  IHEAT  = 1 or  3.  Heat  input  in  (Btu/sec-ft2). 

31-40 

XMINP 

Input  if  IHEAT  = 4.  Mass  flow  rate  of  liquid  impinging 
on  blade  (lbm/sec-ft2). 

41-50 

HIMP 

Input  of  IHEAT  = 4.  Enthalpy  of  liquid  impinging  on 
blade  (Btu/tbm). 

51-60 

DTIME 

Duration  of  ignition  energy  input  in  seconds. 

Input  card 

17  if  IHEAT 

= 2 on  Card  16  input  Nodes  (Card  16)  values. 

Card 

Columns 

Variable 

Definition 

17 

1-2 

INPI(l) 

“I"  location  of  first  node  effected  by  heat  input. 

3-4 

INPJ(l) 

“J”  location  of  first  node  effected  bv  heat  input. 

5-6 

INPI(2) 

7-8 

INPJ(2) 

77-78 

INPK20) 

79-80 

INP.J(20) 

If  more  than  20  points  continue  on  another  card.  All  values  re  right-adjusted  integers. 


I.  OUTPUT  DESCRIPTION 


: 

! 


1.  Initial  external  aerodynamics  at  midspan 


Node 

X-Loc 

Velocity 

Gas 

Press 

Film  Coefficient 


ith  nodal  location 

Chordwise  distance  from  leading  edge,  in. 

Local  air  stream  velocity,  ft/sec 

Air  temperature,  °R 

Air  static  pressure,  PSIA 

Local  heat  transfer  coefficient,  Btu/hr-ft2 — °R 


2.  Graphic  display  of  airfoil  at  each  output  time.  Plot  shows  the  airfoil 
geometry  and  a map  of  the  liquid  metal  region.  Also  included  are  the 
leading  edge  Reynolds  number  and  the  total  aerodynamic  melt  losses  in  the 
axial  and  radial  directions. 


3.  Tabulation  of  melting  velocity  at  each  burning  leading  edge  node  for  each 
output  time. 

X-Loc  X-location  of  node,  in. 

Y-Loc  Y-location  of  node,  in. 

UX  X-component  of  melting  velocity,  ft/sec 

UY  Y-component  of  melting  velocity,  ft/sec 


4.  Tabulation  of  spreading  velocity  at  points  along  the  edge  of  the  liquid  metal 

region  of  each  output  time. 

X-Loc  X-location  of  point,  in. 

Y-Loc  Y-location  of  point,  in. 

UX  X-component  of  melting  velocity,  ft/sec 

UY  Y-component  of  melting  velocity,  ft/sec 


5.  Table  of  temperatures  (°R)  at  each  nodal  location  (i,  j)  at  each  output  time. 
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SECTION  VII 


TEST  CASE 


A representative  test  case  is  enclosed  as  Figures  VII-1  through  VTI-7.  This  case  is  a 
simulation  of  a Ti  8-1-1  airfoil  ( 1.0  in.  chord,  2.0  in.  span,  constant  thickness  of  0.040  in.)  ignited 
at  the  midspan  leading  edge.  The  aerodynamic  environmental  conditions  (temperature,  pressure, 
velocity)  are  1320°R,  100  psia,  and  445  ft/sec.  Figures  VII-1,  -2,  and  -3  are  the  input  data, 
describing  the  geometry,  material  properties,  and  environmental  conditions.  Simulated  output 
data  are  shown  in  Figures  VII-4  through  VII-7,  illustrating  the  output  format. 

Figure  VII-4  shows  the  nodal  temperature  map  in  the  airfoil  0.4  sec  after  ignition.  Burned- 
out  regions  of  the  airfoil  are  shown  in  this  figure  as  being  at  the  melting  temperature.  Figure  VII -5 
shows  the  location  of  the  burning  leading  edge,  the  extent  of  melt  spread  over  the  airfoil,  leading 
edge  regression  rate,  and  the  melt  spreading  velocity.  Figure  VII-6  is  an  unsealed  graphic  display 
printed  on  the  computer  output  page.  Figure  VII-7  is  a Calcomp  plot  of  the  same  data  in  a scaled 
format. 
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Figure  VII-1.  Test  Case  Input  Data 
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-«  f-  ^ r-  r»  ^ r 


e-r^r'-e-r^e-e-e-e-r^i 


8 

00 

s 


H 

Si 

XI 

Z 

CL  ~4 

4 f- 

z 


5 

H- 

5 

£ 

Z 

u 


o 

o 

o 

4* 


u 

z 


O' 

II 


CO 

II 


•o 

II 


\r\ 

ii 


eg  eg  eg  eg'Nrgegegegrgrgegegegegeg^j^jfgrgeg 

■*>  fO  i**!  CO  ^ 1*1  <1  r*l  ^ ■♦I 


eg  <M  rsi  egcgegegcgegrgegrgcgfgcgegcgegcNegfg 


r‘h-fk-rfc^h-r-r*r*r»Kh-r«-h-r»h“r»^h*f'-h* 

<M  eg  eg  eg  «N  eg  N eg  eg  IN  N eg  eg  eg  <N  eg  eg  N >N  *N  eg 

<^."eirOi,41roeO.^.,ei»'r}<'OeOi<li^(*}fOfOf^enf<>(4,>'t) 


fMNNfgN.\Nrg1NNO.M(NNNN,NN,NN(N 
■n  nm^m-n^fn^ro  4 
•O  ■ei  <*)  (*>  <**  <*1  <*» 

H ■— 4 -—4  H H H H H H H fO  H H H ^ *— 4 H H H H H 


fwrK.r'.r-r^f^r-r-r-ooor-r-f^r^f^f^f^r-r^ 

NiNNiMNNNNNOOONNNNMNIMIMN 
.■Oi^f^C".  rO(<k  ro  c^*  U’N  ^ c^  "O  ro  r»%  ro  ^ <■<>  ^ 


NNNNNNNNNOOCNNNN.^rgNNN 

■Ofnr<>f0^c«,>^^^*^»r»r0(nf0*r>  n,'n.-n 

«*>  eo  m >*  m *i  *>  *j  *1  *>  ■n  i*>  *o 


^r'-r-f-r'-Nr-Nf-ooor-f 


^ r*-  r*-  r-  n r-. 


NNNNNNNlNfN/lOlTliMrgNMNNMNN 


•-r^Nr-e-r-r-N^ooor^^-^r-r-e-r-^h* 

.N  N <N  N <N  N r J N N /N  <N  ^ eg  <N  <N  *N  <N  fN  (N  <N  N 

iH^in  > 


<*> 

M 


M 

II 


N 

k r»- 

.7 

.7 

.0 

.0 

.0 

.0 

.0 

.7 

.7 

.7 

.7 

.7 

.7 

.7 

.7 

-N  <N 

N 

•JNIMNNOOI>C 

p ^ eg 

NNN 

*g 

N (N 

^ *> 

.*0 

•o 

^ O fl?> 

a *i  <*\ 

eo  ■n  n 

-0  rO 

"O 

<i  -n  ro 

-0  -o 

^04  4 

t 

-n  m 

“0 

-n  -o 

*-4  r-l 

*-4 

•H  -4 

-4 

-4 

*h  m ^ *1 

«*> 

rO  ^4  -4 

r- 4-4-4 

— i M 

tS 

s 

e-  k 

e- 

r*-  k 

f*- 

h- 

7 

0 

0 

0 

O 

O h-  *» 

7 

7 

7 

7 

7 

7 

3 

• • 

• 

• • 

• 

• 

* * 1 • 

• 

• • • 

• • • 

• 

• • 

eg  eg  eg 

N eg 

N 

eg 

N O ^ ^ 

r- 

O eg  n 

N N N N N N 

E 

•o  •»> 

*\ 

■n  rr\ 

•*\ 

r\  o <r  o 

O' 

6 n <n 

n m *n 

■n 

y%  ir\ 

CL 

^ r\  (*) 

•n  ro 

* *>  M * 4* 

t 

‘jein 

rO  r>  r> 

•*>  -n 

Q, 

-*  .-A 

H 

H 

M 

—4 

*-«  1*0  *e»  *0 

«> 

'O  H H 

*g  -g  -4 

—4 

—4  *-4 

6 

£ 

+m  !*» 

e» 

r-r'-e-e-e-oooo 

o^r^e-e^h-r-^e- 

~a 

(N  eg 

N 

N N 

N 

eg 

NO^'- 

r*- 

© N N 

: n •: 

SI 

N N 

"C3 

'n  >n 

<n 

** 

<^  O 0*  O' 

> 

O^iei 

*\  re>  rev 

•o  *> 

rri 

■n  -n 

2 

•n<^4  t 

t 

m r\  rn 

^ “0  -n 

*0  eo 

^ «-g 

-a 

—4  ^ 4 

~4 

■g^  <»  n 

.*» 

^ -g  — 4 

-<  -*  -4 

-A 

H "A 

-g<N#^^g>  o r*  ® ? o h ^ ^ ® 0>  O -* 

4H-I-IH-I  g j(N*g 


CL 

0c 

£ 


V 1 1 - 5 


TlMfc  = 0.^0000 


j 


z 

o 

£ 

oc 

-J 

< 

t- 

JJ 

X 

c 

o o 

o 

o o 

o 

o 

LU 

o 

r4 

ss 

o 

IA 

c o 

(O  ^ 

0-4 

u 

UJ 

00 

n 

N 

<nj  m 

o 

<✓> 

V 

o 

O O O 

o 

o o 

M 

N 

o 

o 

o © 

o 

o 

o 

h- 

• 

• • • 

• 

• 0 

-J 

H- 

• 

0 

• • 

0 

0 

0 

A. 

o 

o o o 

o 

o o 

a. 

o 

© 

o o 

o 

o 

o 

ID 

JJ 

3 

X 

a 

>- 

*- 

► 

Z 

3 

3 

LL 

aj 

3 

8 

JJ 

UJ 

j? 

j 

■J 

o o o 
o o 5 

o 

o 

8 

5 

jj 

D 

88 

88 

8 

o 

o 

o 

o 

z 

JJ 

NJ  f-  -N 

eg 

jj 

o 

4- 

» © 

CO 

t 

o 

»-4 

./> 

iAO  + 

o 

iTv 

UJ 

A 

4 

h* 

4-  ^ 

4- 

r- 

>*■ 

3 

V 

o 

0-4-4 

»—4 

o o 

r 

N. 

o 

o 

'J  (M 

** 

o 

o 

< 

t- 

• 

0 • • 

0 

• 0 

>- 

• 

• 

• • 

• 

0 

0 

UJ 

LL 

o 

0 

0 

0 

o 

o o 

X 

o 

c 

o o 

o 

o 

o 

-J 

-* 

W 

J 

X 

5 

X 

z 

3 

j 

3 

-4 

< 

•/> 

*- 

CO 

z 

— * 

x 

Z 

o 

o o c 

o 

o o 

o 

Z 

o 

o 

o o 

o 

o 

o 

j 

•“4 

-«  O'  O 

O'  O' 

a. 

»-4 

M 

*4 

O'  o 

»-4 

O' 

4 O'  ~l 

(NJ 

r-  oo 

x> 

4- 

T — 

rsi 

r- 

O0 

0J 

o- 

O'  O'  ^ 

(SI 

vN  N 

H- 

O' 

O' 

O'  M 

N 

N 

(SI 

a 

0 

• • • 

• 

• • 

< 

o 

• 

0 

0 0 

0 

0 

0 

3 

o 

0 

0 

1 

-4 

<-4  -4 

3 

o 

o 

O -J 

<— 1 

r4 

-4 

< 

J 

1 

-J 

1 

>- 

> 

1-^ 

> 

J 

-4 

3 

_> 

-1 

UJ 

-1 

> 

Jj 

> 

Z 

O O M 

o 

■0 

(3 

Z 

o 

o 

o o 

o 

o 

o 

3 

-« 

pi 

S 

z 

H 

o 

00 

8 

*3 

£ 

8 

% 

2 

o » 

» 

> 

a 

o 

o 

oo  -o 

ao 

o 

o 

-4 

J 

o 

O -*  NJ 

-1 

o o 

< 

J 

*■4 

4-  o 

4 

<■> 

*"4 

n 

• • • 

• • 

UJ 

q 

• 

• 

• • 

• 

• 

0 

_j 

□ 

o 

o o o 

o 

o o 

ar 

3 

o 

o 

o o 

O 

o 

o 

JJ 

i 

a. 

i 

r 

X 

A 

X 

Figure  VII -5.  Combustion  Propagation  Data 
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Figure  VII -6.  Approximate  Graphic  Display 
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Figure  VII-7.  Calcomp  Data  Display 
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SECTION  VIII 


< PROGRAM  IDENTIFICATION  AND  REVISION  PROCEDURE 

A.  CCD  NUMBER 

Customer  Computer  Decks  (CCD’s)  are  identified  by  a CCD  number  and  date.  An  example 
is  CCD  1001-0.0,  15  November  1969. 

The  CCD  number  consists  of: 

1.  Basic  number  (first  four  digits) 

2.  Dash  (or  change)  number 

3.  Decimal  (or  addition/correction)  number. 

1.  Basic  Number 

This  four-digit  number  generally  corresponds  to  a given  program.  If  another  method  is 
developed  or  studied  which  does  not  replace  or  supersede  the  original,  a new  four-digit 
number  is  issued. 

2.  Dash  Number 

Major  changes  in  the  program  are  reflected  in  different  dash  numbers.  A change  in 
techniques  or  mathematical  methods  would  produce  a new  dash  number  provided  the  new 
techniques  replace  or  change  the  old  techniques.  For  more  than  nine  changes,  the  dash 
number  shifts  to  letters.  The  dash  number  of  the  original  program  is  zero. 

3.  Decimal  Number 

The  decimal  is  used  for  all  other  program  changes  such  as: 

1.  Adding  new  optional  routines 

2.  A change  to  the  FORTRAN  source  language  due  to  computer  differences 

3.  Correcting  a mistake  in  the  program 

4.  Input  and  output  changes 

5.  Adding  unique  customer  oriented  curves  or  initialization  data. 

In  most  cases  a decimal  number  change  is  made  and  documented  by  writing  an  addendum 
to  the  user’s  manual  without  reprinting  the  user’s  manual.  An  errata  or  addendum  to  the  user’s 
manual  is  used  in  conjunction  with  a Computer  Simulation  Change  Notice  (Section  V of  CCD 
Approval  Form  — Figure  V III  - 1 ) . This  notice  briefly  describes  the  changes  affecting  the  deck  and 
manual. 

Accompanying  the  notice  is  a new  manual  title  page  reflecting  the  new  date  and/or  dash  or 
decimal  changes.  The  SCN  also  includes  the  change  pages  with  black  bars  in  the  right-hand 
margin  opposite  data  changed  and  a new  date  in  the  upper  right-hand  corner. 

Programs  revised  by  the  user  without  written  consent  of  the  supplier  are  the  responsibility 
of  the  user. 
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IV  COMPUTER  SIMULATION  CHANGE  NOTICE  NO 


Program  Title 

Changes  Apply  To: 
(18)  Computer  Deck  No.. 

New  Designations: 
© Computer  Deck  No.. 


.User's  Manual  No.. 


.User's  Manual  No.. 


Type  of  Change:  □ Errata  □ Addendum  Changes  Effect:  □ Deck  □ Manual 


Change: 


Detailed  Reasons  for  Changes  Are  Shown  in  Enclosure.  (~~1 


Approval 

Program  Office.  Engineering  Section  Manager.. 
Engineering  Computer  Applications 


(If  no  change  in  Engineering  Department 
committment,  only  Engineering  Computer 
Applications  sign-off  needed) 


^ Action  Taken  to  Satisfy  Restrictions  - 
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ing Branch  will  sign  m all  cases  For 
non-contract  decks  with  restrictions 
the  Manager  Engineering  Branch 
For  all  other  decks  the  Manager 
Engineering  Branch  may  designate 
his  Section  Manager  to  approve 
transmittal 
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instructions  for  Use  of  Forms  (Continued) 

ill  User  s Manual  -Documentation  Signoff  - 
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documentation  supplied  with  the  deck 
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Figure  VIII- 1.  Computer  Simulation  Change  Notice 
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Date  of  this  notice 
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New  Deck  and  Manual  designations 


APPENDIX  A 


DESCRIPTION  OF  PROGRAM  FUNCTIONS 


Figure  IV- 1 is  a logic  diagram  outlining  the  functional  operation  of  the  program. 
Examination  of  this  diagram  shows  three  major  sections  of  the  program:  (a)  the  input- 
initialization  section,  (b)  the  aerodynamic-heat  transfer  section,  and  (c)  the  combustion-liquid 
metal  section.  In  Section  IV  the  functional  steps  in  the  program  operation  were  outlined.  In  this 
appendix  the  more  complex  functions  will  be  described  in  terms  of  the  equations  and  logic 
incorporated  in  their  operation.  The  functions  are  described  in  the  accompanying  logic  diagrams. 
In  this  appendix,  the  required  equations  are  listed  in  Appendix  B,  and  the  nomenclature  is 
defined  in  Appendix  C. 

The  logic  diagrams  are  prepared  in  symbolic  equation  form.  Thus  ©— 6 means  “Equation 
(1)  from  the  Equation  List  is  used  to  calculate  8,  the  boundary  layer  momentum  thickness.” 
Where  a routine  name,  i.e.,  TEMP,  appears  in  the  logic  flow,  it  refers  to  the  use  of  subroutine 
TEMP  to  yield  defined  information. 

The  aerodynamic  boundary  layer  calculation  will  be  the  first  functional  description.  Prior 
to  this  function  are  the  input  and  initialization  functions  describing  the  environmental  condition, 
initial  temperatures,  properties,  and  initial  geometry  of  the  airfoil. 

A.  BOUNDARY  LAYER  ROUTINE 

Figure  A-l  shows  the  logic  flow.  Entry  is  made  with  the  free  stream  velocity  profile  over  the 
airfoil.  Velocity  may  vary  with  distance,  but  is  considered  symmetrical  on  both  airfoil  surfaces. 
The  function  calculates  the  boundary  layer  growth,  decides  whether  flow  is  laminar  or  turbulent, 
calculates  the  aerodynamic  shear  stress  (CO,  the  convective  heat  transfer  coefficient,  and  the 
mass  transfer  coefficient  (Kxm). 

B.  HEAT  BALANCE  ROUTINE 

As  shown  in  Figure  A-2,  entry  is  made  with  the  metal  node  temperatures,  with  the  surface- 
bonded  oxide  weight  (either  initial  ot  from  prior  calculations),  and  with  the  mass  transfer  and 
heat  transfer  coefficients  from  BOUNDARY  LAYER.  The  surface  oxidation  rate  is  determined 
from  the  selected  kinetic  reaction  model,  based  on  surface  temperature.  The  reaction  heat  input 
is  then  balanced  with  other  modes  of  heat  transfer  to  determine  the  nodal  temperatures 
throughout  the  airfoil.  The  energy  equation  (15)  is  written  in  finite  difference  form  based  on  the 
nodal  arrangement  illustrated  in  Figure  A-3.  The  central  temperatures  Tn,  Tn’  are  defined  as 
those  at  the  beginning  and  end  of  a time  increment.  Equation  (16)  defines  the  time  increment 
parameter  and  Eq.  (17)  is  the  stability  criterion,  defining  the  maximum  allowable  time  increment 
for  stable  calculation. 

After  calculating  the  nodal  temperatures,  they  are  tested  to  determine  whether  melting  has 
occurred  at  any  point  in  the  airfoil.  If  melting  occurs,  the  program  branches  to  the  liquid  metal 
section  to  evaluate  reaction,  flow  and  other  associated  processes.  If  no  melting  occurs,  time  is 
advanced  by  one  increment,  and  the  program  returns  to  perform  another  heat  transfer  analysis. 
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Figure  A-3.  Finite  Difference  Nodal  Arrangement  for  Heat  Transfer  Analysis 


C.  STAGNATION  POINT  ROUTINE 


If  melting  has  occurred,  the  liquid  analysis  section  is  entered  in  the  stagnation  point 
routine,  Figure  A-4,  with  the  airfoil  metal  temperature  thickness,  and  environmental  conditions. 
This  routine  simultaneously  evaluates  the  melting  rate,  oxygen  uptake  rate,  liquid  metal  film 
thickness,  oxygen  concentrations,  and  temperatures  in  the  liquid  metal,  all  at  the  airfoil 
stagnation  edge.  These  calculated  parameters,  illustrated  in  Figure  A-5,  provide  the  initial  data 
for  calculation  of  the  liquid  metal  processes  over  the  airfoil  surface.  An  iterative  procedure  is  used 
to  obtain  the  required  simultaneous  solutions.  On  the  first  pass,  the  liquid  metal  is  assumed  to 
be  at  the  melting  temperature,  allowing  initial  calculation  of  the  various  parameters.  An 
empirical  corelation,  Eq.  ( 19),  provides  an  initial  value  of  the  liquid  film  thickness  at  the  leading 
edge.  Equation  (21)  describes  the  relationship  between  heat  transfer,  reaction,  and  melting  rate, 
and  defines  the  direction  of  surface  regression  at  the  leading  edge.  The  physical  relationship 
between  these  quantities  is  shown  in  Figure  A-6.  In  the  main  iteration  loop,  the  oxygen  absorption 
rate  at  the  liquid  surface  is  defined  by  the  lower  value  of  either  Eq.  (11)  or  (25).  In  the  next  step, 
liquid  metal  temperatures  and  melting  velocity  are  determined  by  subroutine  TEMP,  which  will 
be  separately  described.  On  completion  of  the  iteration,  the  liquid  metal  flow  components  (£,,  h,,) 
and  oxygen  How  rate  in  the  flowing  liquid  (g„)  are  determined.  These  quantities  are  required  as 
initial  data  for  the  following  streamline  analysis.  On  exit  from  the  stagnation  routine,  the  melting 
velocity  (U)  is  tested  to  determine  if  melting  has  ceased. 


D.  TEMPERATURE  ROUTINE 


This  subroutine,  Figure  A-7,  determines  the  temperatures  in  the  liquid  metal  film,  and  the 
melting  velocity  of  the  surface.  The  model  assumes  that  the  liquid  metal  film  is  composed  of  a 
solid-liquid  two-phase  region,  and  a superheated  region,  as  shown  in  Figure  A-8.  The  model  is 
further  developed  on  the  assumption  that  melting  is  distributed  throughout  the  two-phase  region. 
This  subroutine  requires  the  thermodynamic  properties  of  the  alloy,  i.e.,  liquidus  and  solidus 
temperatures  as  functions  of  absorbed  oxygen,  and  enthalpy  as  a function  of  oxygen  content  and 
temperature. 

A-4 


Melt  Flow 


m0  (Oxygen 
Mass  Flow) 


Gas  Flow 


x0  (Oxygen 

Concentration 
in  Gas) 


Liquid 

Solid 


. (Film  Thickness) 

■•r: y-  r . * rv.  ••  \\ ", 

U (Melting  Velocity) 


(Oxygen  : 

Concentration  in  Melt) 


Leading 
' Edge 
Thickness 
(A)  \ 


Figure  A-5.  Stagnation  Point  Region  of  Burning  Airfoil 
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Figure  A -6.  Heat  Flow  and  Velocity  Relationship  at  Burning  Leading  Edge 


AD-A072  529 


UNCLASSIFIED 


PRATT  ANO  WHITNEY  AIRCRAFT  6R0UP  WEST  PALM  BEACH  FL  6—ETC  F/6  11/6 
TITANIUM  CMBUSTION  RESEARCH  PROGRAM  AND  USER'S  MANUAL  FOR  DECK  —ETC(U) 
MAR  79  MR  GLICKSTEIN  F33615-76-C-5041 


PWA-FR-10748 


AFML  -TR-79-4001 


2 OF  2 
*8^2529 


9-79 


! 


Figure  A -8.  Local  Velocity  Profiles  and  Nomenclature  Used  in  Melt  Flow 
Analysis 

E.  STREAMLINE  ROUTINE 

On  completion  of  evaluation  of  the  leading  edge  stagnation  points,  the  program  performs  an 
adjustment  of  the  time  increment.  The  recession  of  the  leading  edge  due  to  melting  during  one 
time  increment  is  calculated  based  on  the  time  increment  defined  by  stability  requirements  (Eq. 
17)  and  the  melting  velocity  calculated  at  the  leading  edge.  If  the  recession  distance  exceeds  the 
node  spacing,  the  time  increment  is  reduced  to  make  the  two  distances  approximately  equal. 

Entry  is  then  made  into  the  streamline  routine.  Figure  A-9,  for  each  node  point  at  the 
leading  edge,  and  the  liquid  flow  processes  are  evaluated  along  the  streamline  paths  starting  at 
each  leading  edge  node,  as  shown  in  Figure  IV-2.  The  differential  liquid  metal  transport  equations 
(Eqs.  35,  36,  37)  are  integrated  with  an  implicit-explicit  Taylor  series  to  determine  the  streamline 
growth  over  the  surface.  An  iterative  procedure  is  then  employed  to  evaluate  the  remaining  liquid 
parameters  along  the  streamlines  (oxygen  content,  film  thickness,  temperature,  etc.).  During  this 
process  a correction  is  made  to  the  calculated  liquid  film  thickness  to  account  for  the  effects  of 
streamline  interaction  (crowding  together  .>r  spreading  apart).  Calculation  along  each  streamline 
continues  until  the  downstream  edge  of  the  liquid  region  is  reached. 
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F.  FLOWING  MELT 


The  advance  of  molten  liquid  over  the  airfoil  surface  is  governed  by  a balance  of  forces; 
aerodynamic  drag  and  centrifugal  forces  driving  the  liquid,  and  surface  tension  retarding  its 
motion,  as  shown  in  Figure  A-10.  As  the  liquid  approaches  the  end  of  the  streamline,  its  velocity 
slows  and  the  stream  thickens.  Balancing  forces  at  the  end  of  the  stream  yields  Eq.  (44),  which 
is  combined  with  Eqs.  (30)  and  (45)  to  yield  the  spreading  velocity  (Ve),  the  rate  at  which  the 
burning  molten  metal  spreads  across  the  airfoil. 

At  the  completion  of  the  streamline  calculations,  the  program  sums  the  total  loss  of  molten 
metal  to  the  airstream.  The  change  in  airfoil  geometry  is  calculated  from  the  burning  velocities 
(U)  over  the  airfoil,  and  the  program  returns  to  the  aerodynamic  section  to  re-evaluate  the  airfoil 
temperatures. 
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EQUATION  LIST 


u»*’  - °n4Jf  U*  dx 

V.  Ui 

(Laminar  Flow) 

(1) 

e ( Y=  u-‘  [c,  + 0.16  £ UJ  dx  ] 

(Turbulent  Flow) 

(2) 

Re0  = K = — ^ 

v,  v„  dx 

(3) 

„ func  (K) 

L’  Red 

(Laminar) 

(4a) 

C,  = 0.0256  Re0'“ 

(Turbulent) 

(4b) 

N o PrR ed  V Cf/2 

“ 5Pr  + 5 In  (5Pr  + 1)  + n/  2/Cf  - 14 

(Laminar) 

(5a) 

Nu* = nr 

(5b) 

Nufl  = *«  Pf  * Cf 

(Turbulent) 

(6) 

Nu„  = 1.14  Re?  Pr0,7‘ 

(7a) 

Nuo 

(7b) 

Re0=<^ 

(7c) 

f(>)  = [ 0.0242  e'2,>'’r/2)  + cos  y ]°®  + 0.05  sin  y 

(8a) 

0 < 7 < ir/2 

f(>)  = |[  0.0242  e •2(t  - *m  + cos  7 ]°“  + 0.05  sin  7 } ( 3.5  -»•«»*  ) 

(8b) 

t/2  < 7 < w 

Nu  = NUc.lc  [l  - (l  -f(y))e  ^] 

(9) 

arg  = -0.016  Red  [If  arg  < -180]  arg  = -180 


Kxm  = gP 


(10) 


B-l 


= 7.17  A'0*  + 0.00673A 

pU(0,  - $) 

Mi  (T,  4>)  = vm  p [ y j-j‘  ] , T,  < T < T, 

T,  = j (T,  + T,) 
f„  = 0 

ho  = 0 

F = V P + h1  , 

B 


df  I7  , 

fa  = U cos  i - £f 

(35) 

, t8 1 

6 Mi 

(35a) 

-£■  = U sin  0 - f h 

(36) 

. 08*  t0  8 1 

12pi  6mi 

(36a) 

^ = m0  - £0,  pF 

(37) 

f(s,k,i)  - f (s  - l.i)  + ^f-  [f.  (s  - l,i  - 1)  + f„(s  - l,i)] 

(38) 

h (s,k,i)  = h (s  - l.i)  +4p  [h,  (s  - l.i  - 1)  + h,(s  - l.i)] 

(39) 

g (s,k,i)  = g (s  - l.i)  +—  [g,  (a  - l.i  - 1)  + g,(s  - l.i)] 

g (s,k,i)  = g (a  - l.i)  + 4r  [g.  (a  - l.i  - 1)  + g.(a  - l.i)] 

<U  s <p  + M TT 

(40) 

(41a) 

W (2f»  + 3hJ  + hf  y-  ) 

(41b) 

-uJu 

v'-  (M\  81 

(41c) 

V 2m,  / 

0 = pRoj* 

(41d) 

Aicorr  = [ A6  ]w- const  + 'jjj"  * AS 

(42) 

*-£ -r-^ 

(43) 

(44) 

V . Z 
* \ 

(45) 
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APPENDIX  C 


LIST  OF  SYMBOLS 


Notation 


A|y 

CD 

CF 

D 

Don 

F 

Fr 

f 


g 

h 

he 

H 

H„, 

H„ 

?J 

in 


k 

K 

K,m 


j+  i- 
1 Xt  1 X 


m0 


n 

A*.  q* 

R 

ft 

S.J.I 


8 


t 

T 

T. 


= surface  area  of  nodal  mass 
= aerodynamic  drag  coefficient 
= coefficient  of  friction 
= oxygen  diffusivity  in  metal 
= binary  oxygen -nitrogen  diffusion  coefficient 
= melt  flow  scalar 
= radiation  configuration  factor 
= chordwise  component  of  melt  flow 
= oxygen  convective  flow  in  melt 
= radial  component  of  melt  flow 
= convective  heat  transfer  coefficient 
= enthalpy  of  metal 

= heat  of  reaction  per  unit  mass  of  oxygen 
= latent  heat  of  fusion 

= unit  vectors  in  chordwise  and  radial  directions 
= unit  vector  normal  to  surface 
= thermal  conductivity 
= boundary  layer  shape  factor 
= mass  transfer  coefficient 
= nodal  x-spacing 
= nodal  z-spacing 
= mass  of  oxygen  bonded  on  surface 
= oxygen  mass  flux 
= direction  normal  to  surface 
= surface  heat  flux  vector,  scalar 
= internal  conduction  vector,  scalar 
= leading  edge  aspect  ratio 
= radius  to  mid-span  of  blade 
= streamline  position,  location  index,  and  time  index 
= distance  along  streamline 
= time 

= temperature  in  melted  metal 
= gas  temperature 
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LIST  OF  8YMBOL8  (CONTINUED) 


Notation 


TS,  T; 
T+„  T-, 
U 

u„  u, 
u. 

u,  v 
V,  V 
Ve 
w 

X,Y,Z 

Xe,Ze 

xt 

a 

0 

y 

r 

6 

do 

6, 

A 

A0 

r 

e 

0 

X 

A 

K 

M 

V 

Vo 

Vm 


= melting  temperature 
= node  temperature 
= temperatures  of  x-adjacent  nodes 
= temperature  of  z-adjacent  nodes 
= melting  velocity  of  surface 
= x and  z components  of  melting  velocity 
= free-stream  air  velocity 

= chordwise  and  radial  velocity  components  (melt) 
= liquid  metal  velocity  vector,  scalar 
= spreading  velocity  of  liquid  metal 
= width  of  streamline 

= chordwise,  normal,  radial  coordinate  system 
= end  of  liquid  metal  region 
= laminar-turbulent  transition  point 
= fraction  of  nodal  surface  covered  by  melt 
= centrifugal  body  force  per  volume 
= leading  edge  yaw  angle 
= mass  transport  ratio 
= liquid  metal  film  thickness 
= thickness  of  liquid  at  leading  edge 
= thickness  of  two-phase  region 
= local  thickness  of  airfoil 
= leading  edge  thickness 

= component  of  surface  tension  parallel  to  surface 
= boundary  layer  momentum  thickness 
= Fourier  modulus 

= thickness  of  liquid  metal  leading  edge 
= mass  transfer  Peclet  number 
= thermal  diffusivity  of  solid 
= dynamic  viscosity  of  liquid  metal 
= kinematic  viscosity  of  liquid  metal 
= kinematic  viscosity  of  air 
= kinematic  viscosity  at  liquidus  temperature 
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LIST  OF  SYMBOLS  (Concluded) 


Notation 

i 

p 

p» 

<r 

T 

To 

</> 

X0»o 

^o»oo 

u 

NU9 

Pr 

Reg 

Sc 

St 

Subscripts 


= fraction  of  melt  flow  loss  jom  surface  per  unit  length 
= metal  density 
= air  density 

= Stefan-Boltzmann  constant 
= aerodynamic  shear  stress  in  chordwise  direction 
= aerodynamic  shear  stress  in  radial  direction 
= local  oxygen  mass-fraction  in  liquid  metal 
= bulk  (mean)  oxygen  mass-fraction  in  liquid  metal 
= oxygen  mol-fraction  in  air  at  surface 
= oxygen  mol-fraction  in  air  free-stream 
= streamline  angle  relative  to  chord 
= rotational  speed  of  blade 
= Nusselt  number  for  heat  transfer 
= Prandtl  number 

= Reynolds  number  based  on  momentum  thickness 
- Schmidt  number 
= Stanton  number  for  mass  transport 


= interface  of  liquid  and  two-phase  regions 
= solid 
= liquid 

= liquid-solid  interface 
= liquid  metal-air  interface 
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